Search NASA⌕ Search

SEARCH · Search NASA

Results for “HORIZONTAL TEMPERATURE DISTRIBUTIONS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

210 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

TMI Rain Rate Estimation Over Land and Ocean Utilizing Convective and Stratiform Discrimination

Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) radiometer brightness temperature data in the 85 GHz channel (T85) reveal distinct local minima in a regional map containing a Mesoscale Convective System (MCS). This is because of relatively small footprint size (approximately 5.5 km) and strong extinction properties in this channel of the TMI. A map of rain rate for that region, deduced from simultaneous measurements made by the Precipitation Radar (PR) on board the TRMM satellite, reveals that these T85 minima, produced by scattering, correspond to local PR rain maxima. Utilizing the PR rain rate map as a guide, we infer from TMI data the presence of three different kinds of thunderstorms or Cbs. They are young, mature, and decaying Cbs that have a scale of about 20 km on the average. Two parameters enable us to infer these three kinds of Cbs objectively: a) the magnitude of scattering depression deduced from local T85 minima and b) the mean horizontal gradient of T85 around such minima. Knowing the category of a given Cb, we can estimate the rain rate associated with it. Such estimation is done with the help of relationships linking T85 minimum to rain rate in each Cb type. Similarly, a weak background rain rate in all the areas where T85 is less than 260 K is deduced with another relationship linking T85 to rain rate. In our rain retrieval model, this background rain constitutes the stratiform rain where the Cbs are absent. Initially, these relationships are optimized or tuned utilizing the PR and TMI data of a few MCS events. After such tuning, the model is applied to independent MCS cases. The areal distribution of light (1-10 mm/hr), moderate (10-20 mm/hr), and intense (> 20 mm/hr) rain rates are retrieved satisfactorally. Accuracy in the estimates of the light, moderate and intense rain areas and the mean rain rates associated with such areas in these independent MCS cases is on the average about 15%. Taking advantage of this ability of our retrieval method, one could derive the latent heat input into the atmosphere over the 760 km wide swath of the TMI radiometer in the tropics.

Prabhakara, C.↗
NASA NTRS1993

Tidal and atmospheric forcing of the upper ocean in the Gulf of California. 2: Surface heat flux

Satellite infrared imagery and coastal meteorological data for March 1984 through February 1985 are used to estimate the net annual surface heat flux for the northern Gulf of California. The average annual surface heat flux for the area north of Guaymas and Santa Rosalia is estimated to be +74 W/sq m for the 1984-1985 time period. This is comparable to the +20-50 W/sq m previously obtained from heat and freshwater transport estimates made with hydrographic surveys from different years and months. The spatial distribution of the net surface heat flux shows a net gain of heat over the whole northern gulf. Except for a local maximum near San Esteban Island, the largest heat gain (+110-120 W/sq m) occurs in the Ballenas and Salsipuedes channels, where strong tidal mixing produces anomalously cold sea surface temperatures (SSTs) over much of the year. The lowest heat gain occurs in the Guayamas Basin (+40-50 W/sq m), where SSTs are consistently warmer. In the relatively shallow northern basin the net surface heat flux is farily uniform, with a net annual gain of approxmately +70 W/sq m. A local minimum in heat gain (approximately +60 W/sq m) is observed over the shelf in the northwest, where spring and summer surface temperatures are particularly high. A similar minimum in heat gain over the shelf was observed in a separate study in which historical SSTs and 7 years (1979-1986) of meteorological data from Puerto Penasco were used to estimate the net surface heat flux for the northern basin. In that study, however, the heat fluxes were higher, with a gain of +100 W/sq m over the shelf and +114 W/sq m in the northern basin. These larger values are directly attributable to the higher humidities in the 1979-1986 study compared to the 1984-1985 satellite study. High humidities reduce evaporation and the associated latent heat loss, promoting a net annual heat gain. In the norther Gulf of California, however, tidal mixing appears to play a key role in the observed gain of heat. Deep mixing in the island region produces a persistent pool of cold water which is mixed horizontally by the large-scale circulation, lowering surface temperatures over most of the northern gulf. These cold SSTs decrease evaporation by reducing the saturation vapor pressure of the overlying air. As a result, heat loss is substantially reduced, even when humidities are low. By removing heat from the surface, tidal mixing alters the time scale of air-sea interaction and reduces or possibly even inhibits the formation of deep water masses via convection.

Paden, Cynthia A.↗
NASA NTRS1999

Environmental Forcing of Super Typhoon Paka's (1997) Latent Heat Structure

The distribution and intensity of total (i.e., combined stratified and convective processes) rainrate/latent heat release (LHR) were derived for tropical cyclone Paka during the period 9-21 December, 1997 from the F-10, F-11, F-13, and F-14 Defense Meteorological Satellite Special Sensor Microwave/Imager and the Tropical Rain Measurement Mission Microwave Imager observations. These observations were frequent enough to capture three episodes of inner core convective bursts that preceded periods of rapid intensification and a convective rainband (CRB) cycle. During these periods of convective bursts, satellite sensors revealed that the rainrates/LHR: 1) increased within the inner eye wall region; 2) were mainly convectively generated (nearly a 65% contribution), 3) propagated inwards; 4) extended upwards within the middle and upper-troposphere, and 5) became electrically charged. These factors may have caused the eye wall region to become more buoyant within the middle and upper-troposphere, creating greater cyclonic angular momentum, and, thereby, warming the center and intensifying the system. Radiosonde measurements from Kwajalein Atoll and Guam, sea surface temperature observations, and the European Center for Medium Range Forecast analyses were used to examine the necessary and sufficient condition for initiating and maintaining these inner core convective bursts. For example, the necessary conditions such as the atmospheric thermodynamics (i.e., cold tropopause temperatures, moist troposphere, and warm SSTs [greater than 26 deg]) suggested that the atmosphere was ideal for Paka's maximum potential intensity (MPI) to approach super-typhoon strength. Further, Paka encountered weak vertical wind shear (less than 15 m/s ) before interacting with the westerlies on 21 December. The sufficient conditions, on the other hand, appeared to have some influence on Paka's convective burst, but the horizontal moisture flux convergence values in the outer core were weaker than some of the previously examined tropical cyclones. Also, the upper tropospheric outflow generation of eddy relative angular momentum flux convergence was 4D much less than that found during moderate tropical cyclone/trough interaction. These results indicated how important the external necessary condition and the internal forcing (i.e., CRB cycle) were in generating Paka's convective bursts as compared to the external sufficient forcing mechanisms found in higher latitude tropical cyclones. Later, as Paka began to interact with the westerlies, both the necessary (i.e., strong vertical shear and colder SSTs) and sufficient (i.e., dry air intrusion) external forcing mechanisms helped to decrease Paka's rainrate.

Rodgers, Edward B.↗
NASA NTRS1999

A TRMM Microwave Radiometer Rain Rate Estimation Method with Convective and Stratiform Discrimination

Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) radiometer brightness temperature data in the 85 GHz channel (T85) reveal distinct local minima (T85min) in a regional map containing a Mesoscale Convective System (MCS). This is because of relatively small footprint size (approximately 5.5 km) and strong extinction properties in this channel of the TMI. A map of surface rain rate for that region, deduced from simultaneous measurements made by the Precipitation Radar (PR) on board the TRMM satellite, reveals that these T85(sub min), produced by scattering, correspond to local PR rain maxima. Utilizing the PR rain rate map as a guide, we infer empirically from TMI data the presence of three different kinds of thunderstorms or Cbs. These Cbs are classified as young, mature, and decaying types, and are assumed to have a scale of about 20 km on the average. Two parameters are used to classify these three kinds of Cbs based on the T85 data: a) the magnitude of scattering depression deduced from local T85(sub min) and b) the mean horizontal gradient of T85 around such minima. Knowing the category of a given Cb, we can estimate the rain rate associated with it. Such estimation is done with the help of relationships linking T85min to rain rate in each Cb type. Similarly, a weak background rain rate in all the areas where T85 is less than 260 K is deduced with another relationship linking T85 to rain rate. In our rain retrieval model, this background rain constitutes stratiform rain where the Cbs are absent. Initially, these relationships are optimized or tuned utilizing the PR and TMI data of a few MCS events. After such tuning, the model is applied to independent MCS cases. The areal distribution of light (1-10 mm/hr), moderate (10-20 mm/hr), and intense (>= 20 mm/hr) rain rates are retrieved satisfactorily. Accuracy in the estimates of the light, moderate, and intense rain areas and the mean rain rates associated with such areas in these independent MCS cases is on the average about 15 %. Taking advantage of this ability of our retrieval method, one could derive the latent heat input into the atmosphere over the 760 km wide swath of the TMI radiometer in the tropics.

Prabhakara, C.↗
NASA NTRS1994

Equatorial wave activity derived from fluctuations in observed convection

The spectrum of equatorial wave activity propagating vertically into the stratosphere is calculated from high-resolution imagery of the global convective pattern. Synoptic Global Cloud Imagery (GCI), constructed from six satellites simultaneously observing the earth, is used to diabatically force the linearized primitive equations. Having resolution of 0.5 deg and 3 h, that imagery captures the dominant scales of organized convection, including several harmonics of the diurnal cycle. Its global coverage with high space-time resolution allows the GCI to represent heating variability and dynamical behavior excited by it over a wide range of scales. The dynamical response above the heating is evaluated globally in terms of a space-time spectrum of Hough modes, one which includes planetary-scale Kelvin waves, Rossby waves, and gravity waves down to the resolution of the GCI. The geopotential response, which is indicative of temperature fluctuations observed by satellite, is very red in frequency. Therefore, planetary-scale waves with periods longer than two days dominate the spectrum of geopotential, while high-frequency gravity waves make a comparatively small contribution. Some 80% of the geopotential variance is accounted for by the Kelvin and gravest-symmetric Rossby modes, while the Rossby-gravity mode is comparatively weak. In horizontal eddy motion, the excited wave spectrum is still dominated by planetary-scale components. However, meridional wind fluctuations associated with the Rossby-gravity mode have variance comparable to that of zonal wind fluctuations associated with the Kelvin mode, even though the Rossby-gravity mode is nearly invisible in the geopotential response. Estimates of tropospheric heating lead to amplitudes and propagation characteristics that are broadly consistent with satellite and radiosonde observations of wave activity in the lower stratosphere. The space-time spectrum of EP flux is significantly whiter than the response in either geopotential or motion. Gravity waves of small scale and high frequency carry a large fraction of the upward flux. Although it dominates eastward variance of geopotential and motion, the Kelvin mode carries only about 50% of the eastward EP flux at phase speeds of 20-40 m/s and only 35% of the total eastward flux transmitted to the stratosphere. The remainder is carried by the gravity wave spectrum, which carries nearly all of the westward flux at phase speeds greater than 20 m/s. The gravity wave spectrum also contributes significantly at phase speeds of 10-20 m/s, where only 25% of the flux is accounted for by zonal wavenumbers less than 20. The broad nature of the gravity wave spectrum suggests its absorption at critical levels will be distributed over a deep layer of the middle atmosphere.

Bergman, John W.↗
NASA NTRS2002

A New Formation Mechanism for the Hottest Horizontal-Branch Stars

Hot subluminous stars lying up to 0.7 mag below the extreme horizontal branch (EHB) are found in the ultraviolet color-magnitude diagrams (CMDs) of both omega Cen and NGC 2808. In order to investigate the origin of these subluminous stars, we have constructed a detailed set of evolutionary sequences that follow the evolution of low-mass stars continuously from the zero-age main sequence through the helium-core flash to the HB for a wide range in the mass loss along the red-giant branch (RGB). Stars with the largest mass loss evolve off the RGB to high effective temperatures before igniting helium in their cores. Our results indicate that the subluminous EHB stars, as well as the high temperature gap along the EHB of NGC 2808, can be explained if these stars undergo a late helium-core flash while descending the white-dwarf cooling curve. Under these conditions the convection zone produced by the main helium flash will penetrate into the stellar envelope, thereby mixing most, if not all, of the envelope hydrogen into the hot helium-burning interior, where it is rapidly consumed. This phenomenon is analogous to the 'born-again' scenario for producing hydrogen-deficient stars following a very late helium-shell flash. This 'flash mixing' of the envelope during a late helium-core flash greatly enhances the envelope helium and carbon abundances and, as a result, leads to a discontinuous increase in the HB effective temperature. We argue that the hot HB gap observed in NGC 2808 is associated with this theoretically predicted dichotomy in the RB properties. Using new helium- and carbon-rich stellar atmospheres, we show that the changes in the envelope abundances due to flash mixing will suppress the ultraviolet flux in the spectra of hot EHB stars. We suggest that such changes in the emergent spectral energy distribution are primarily responsible for explaining the hot subluminous EHB stars in omega Cen and NGC 2808. Moreover, we demonstrate that models without flash mixing lie, at most, only approximately 0.1 mag below the EHB, and hence fail to explain the observations.

Sweigart, Allen V.↗
NASA NTRS2012

Projections of Declining Surface-Water Availability for the Southwestern United States

16 of the CMIP5 models had all the data needed for this work for at least one simulation that was continuous from 1950 to 2040. Details of the models analyzed here are provided in Table S1. The model data analyzed here are available at http://strega.ldeo.columbia.edu:81/expert/home/.naomi/.AR5/.v2/.historical:rcp85/.mmm16/ a. Assessing the climatology of the models Despite increases in horizontal resolution of many models compared to their CMIP3 counterparts none of these models can adequately resolve the topography of the south west United States, such as the Sierra Nevada and Rocky Mountains and the associated orographic precipitation. This requires that caution be used when interpreting the results presented here. To assess the ability of the models to simulate the current hydroclimate, in Figure S1 we show the observed (from the Global Precipitation Climatology Centre gridded rain gauge data, (1)) monthly climatology of precipitation and the same for all the models and the multimodel mean for the California-Nevada, Colorado headwaters and Texas regions. The GPCC data uses rain gauges only and interpolates to regular grids of which we used the 1◦ by 1◦ one. Details of the data set can be found in (2). While the models apparently overestimate precipitation in California and Nevada the seasonal cycle with wet winters and dry summers is very well represented. It is also possible that the rain gauge observations are biased low by inadequately sampling the higher mountain regions. How ever the models might also be expected to underestimate orographic precipitation due to inadequate horizontal resolution. The 25 models are also too wet in the Colorado headwaters region but correctly represent the quite even distribution though the year. The bimodal distribution of precipitation in Texas, with peaks in May and September, and the absolute amounts, are well modeled but with the September peak too weak. The positive precipitation bias translates into a positive runoff bias for the Colorado headwaters as also shown in Figure S1. Here the observed runoff values are taken from simulations of the Variable Infiltration Capacity (VIC) land surface-hydrology model (3) forced by observed meteorology (5) that were conducted as part of the North American Land Data Assimilation System project phase 2 ( (NLDAS-2), http://www.emc.ncep.noaa.gov/mmb/nldas/. Runoff for California-Nevada is better simulated but there is a positive bias over Texas despite no strong precipitation bias. To check whether regional climate models better simulate P and runoff in these regions we analyzed the historical simulation with the Regional Climate Model version 3 driven by the National Centers for Environmental Prediction-Department of Energy Reanalysis 2 available from the North American Regional Climate Change Assessment Program (http://www.narccap.ucar.edu). This model configuration retained these biases in P and runoff although they were reduced in amplitude. Given these varying biases we plot P and P − E changes in actual values but apply the simplest bias correction possible to the runoff and soil moisture values and show the modeled changes in terms of percentages of the 20th Century model climatologies. A thorough assessment of the simulation of North American climate in CMIP5 models is conducted in Sheffield at al. (North American Climate in CMIP5 Experiments. Part I: Evaluation of 20th Century Continental and Regional Climatology, manuscript submit ted to J. Climate, available at http://www.climate.noaa.gov/index.jsp?pg=./cpo pa/ mapp/cmip5 publications.html). Sheffield et al. analyze the climatology of precipitation, surface air temperature, low level winds, moisture fluxes, runoff etc. and conclude that the main features of the hydrological cycle, including characteristics of the atmospheric moisture balance and its seasonality, are captured in the CMP5 models subject to biases in total precipitation amounts. We chose to use all available models instead of selecting some and rejecting others based on an assessment of model realism. This is in accord with the suggestions of Mote et al. for CMIP3 (4) but future work needs to revisit this matter for the case of the CMIP5 ensemble.

Seager, Richard↗
NASA NTRS2001

The Origin of Hot Subluminous Horizontal-Branch Stars in Omega Centauri and NGC 2808

Hot subluminous stars lying up to 0.7 mag below the extreme horizontal branch (EHB) are found in the ultraviolet color-magnitude diagrams of both omega Cen (DCruz et.al. 2000) and NGC 2808 (Brown et al. 2001). In order to investigate the origin of these subluminous stars, we have constructed a detailed set of evolutionary sequences that follow the evolution of low-mass stars continuously from the zero-age main sequence through the helium-core flash to the HB for a wide range in the mass loss along the red-giant branch (RGB). Stars with the largest mass loss evolve off the RGB to high effective temperatures before igniting helium in their cores (Castellani & Castellani 1993). Our results indicate that the location of the subluminous EHB stars, as well as the high temperature gap along the EHB of NGC 2808, can be explained if these stars undergo a late helium-core flash while descending the white-dwarf cooling curve. Under these conditions the convection zone produced by the main helium flash will penetrate into the stellar envelope, thereby mixing most, if not all, of the envelope hydrogen into the hot helium-burning interior, where it is rapidly consumed (Sweigart 1997). This phenomenon is analogous to the "born-again" scenario for producing hydrogen-deficient stars following a very late helium-shell flash. This "flash mixing" of the envelope during a late helium-core flash greatly enhances the envelope helium and carbon abundances and, as a result, leads to a discontinuous increase in the HB effective temperature. We argue that the hot HB gap observed in NGC 2808 is associated with this theoretically predicted dichotomy in the HB properties. Using new helium- and carbon-rich stellar atmospheres, we show that the changes in the envelope abundances due to flash mixing will suppress the ultraviolet flux in the spectra of hot EHB stars. We suggest that such changes in the emergent spectral energy distribution are primarily responsible for explaining the hot subluminous EHB stars in omega Cen and NGC 2808. Moreover, we demonstrate that models without flash mixing lie, at most, only approximately 0.1 mag below the EHB, and hence fail to explain the observations.

Sweigart, Allen V.↗
NASA NTRS1999

A Study of Direct and Cloud-Mediated Radiative Forcing of Climate Due to Aerosols

The Intergovernmental Panel on Climate Change (IPCC) has reported that in the southeastern US and eastern China, the general greenhouse warming due to anthropogenic gaseous emissions is dominated by the cooling effect of anthropogenic aerosols. To verify this model prediction in eastern China and southeastern US, we analyzed regional patterns of climate changes at 72 stations in eastern China during 1951- 94 (44 years), and at 52 stations in the southeastern US during 1949-94 (46 years) to detect the fingerprint of aerosol radiative forcing. It was found that the mean rates of change of annual mean daily, maximum, minimum temperatures and diurnal temperature range (DTR) in eastern China were 0.8, -0.2, 1.8, and -2.0 C/100 years respectively, while the mean rates of change of annual mean daily, maximum, minimum temperatures and DTR in the southeastern US were -0.2, -0.6, 0.2, and -0.8 C/100 years, respectively. This indicates that the high rate of increase in annual mean minimum temperature in eastern China results in a slightly warming trend of daily temperature, while the high rate of decrease in annual mean maximum temperature and low rate of increase in annual mean minimum temperature lead to the cooling trend of daily temperature in the southeastern US. We found that the warming from the longwave forcing due to both greenhouse gases and aerosols was completely counteracted by the shortwave aerosol forcing in the southeastern US in the past 46 years. A slightly overall warming trend in eastern China is evident; winters have become milder. This finding is explained by hypothesizing that increasing energy usage during the past 44 years has resulted in more coal and biomass burning, thus increasing the emission of absorbing soot and organic aerosols in eastern China. Such emissions, in addition to well-known Asia dust and greenhouse gases, may be responsible for the winter warming trend in eastern China that we have reported here. The sensitivity of aerosol radiative properties to aerosol composition, size distribution, relative humidity (RH) is examined for the following aerosol systems: inorganic and organic ions (Cl-, Br-, NO3 -, SO4 2-, Na+, NH4 +, K+, Ca2+, Mg2+, HCOO-, CH3COO-, CH3CH2COO-, CH3COCOO-, OOCCOO2-, MSA-1); water-insoluble inorganic and organic compounds (elemental carbon, n-alkanes, SiO2, Al2O3, Fe2O3 and other organic compounds). The partial molar refraction method was used to calculate the real part of the refractive index. It was found that the asymmetry factor increased by approximately 48% with the real part varying from 1.40 to 1.65, and the single scattering albedo decreased by 24% with the imaginary part varying from -0.005 to -0.1. The asymmetry factor increased by 5.4 times with the geometric standard deviation varying from 1.2 to 3.0. The radiation transmission is very sensitive to the change in size distribution; other factors are not as significant. To determine the aerosol direct radiative forcing (ADRF), the aerosol optical depth (AOD) values at the three operational wavelengths (415, 500 and 673 nm) were determined at a regionally representative site, namely, Mt. Gibbs (35.78 deg N, 82.29 deg W, elevation 2006 m) in Mt. Mitchell State Park, NC, and a site located in an adjacent valley (Black Mountain, 35.66 deg N, 82.38 deg W, elevation 951 m) in the southeastern US. The two sites are separated horizontally by 10 km and vertically by 1 km. It was found that the representative total AOD values at 500 nm at the valley site for highly polluted (HP), marine (M) and continental (C) air masses were 0.68 +/- 0.33, 0.29 +/- 0.19 and 0.10 +/- 0.04, respectively. A search-graph method was used to retrieve the columnar size distribution (number concentration N, effective radius reff and geometric standard deviation=?g) from the optical depth observations at three operational wavelengths. The ground albedo, single scattering albedo and imaginary part of the refractive index were calculated using a mathematically unique procedure involving a Mie code and a radiative transfer code in conjunction with the retrieved aerosol size distribution, AOD, and diffuse-direct irradiance ratio. It was found that N, r(eff) and sigma(g) were in the ranges of 10 to 1.7 x 10(exp 4)/cubic cm, 0.09 to 0.68 micrometers and 1.12 to 2.95, respectively. The asymmetry factor and single scattering albedo were in the ranges of 0.63 to 0.75 and 0.74 to 0.97 respectively. The ground albedo for the forested terrain and imaginary part of refractive index were found to be in the ranges of 0.06 to 0.29 and 0.005 to 0.051 respectively. On the basis of these aerosol radiative properties obtained at the research sites and computations using the Column Radiation Model (CRM) of National Center of Atmospheric Research (NCAR) Community Climate Model (CCM3), it was found that the average cloud-free 24-hour ADRF values were -13 +/- 8, -8 +/- 3, -33 +/- 16 W/square m for marine, continental, and polluted air masses, respectively. On the assumption that the fractional coverage of clouds is 0.61, it was estimated that the annual mean ADRF was 7 +/- 2 W/square m in the southeastern US. The review with respect to the current knowledge of organic acids shows that aerosol formate and acetate concentrations range from 0.02 to 5.3 nmol/cubic m and from 0.03 to 12.4 nmol/cubic m respectively, and that between 34% to 77% of formate and between 21% to 66% of acetate are present in the fine fraction of aerosols. It was found that although most (98-99%) of these volatile organic acids were present in the gas phase, their concentrations in the aerosol particles were sufficient to make them a good candidate for cloud condensation nuclei (CCN). It is hypothesized that organic acids are at least one of the primary sources of CCN in the atmosphere due to their ubiquitous presence in the troposphere, especially over the continental forested areas. The results of our measurements at Palmer Station, Antarctica show that the daily average CCN concentrations at 0.3% and 1% supersaturations ranged from 0.3 to 160/cubic cm and from 4 to 168/cubic cm, respectively, during the period from 17 January to 26 February, 1994. New evidence for substantial and definitive CCN enhancement near and within cloud has been observed at Mt. Mitchell, North Carolina. The results show that the average monthly CCN concentrations were 460 +/- 217, 386 +/- 286, 429 +/- 228 and 238 +/- 134/cubic cm for in-cloud, overcast, clear and rainy conditions, respectively. The typical CCN spectra show that there were a lot of small CCN produced and the ion concentrations (especially H+ and SO4 2-) were very high during the CCN enhancement period. The significantly positive correlation between black carbon (BC) and CCN at 1% supersaturation indicates that a percentage of the BC measured at the site may be in the form of an internal mixture and participated in the formation of CCN.

Yu, Shao-Cai↗
NASA NTRS1961

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗
NASA NTRS2016

Photospheric Current Spikes and Their Possible Association with Flares - Results from an HMI Data Driven Model

A data driven, near photospheric magnetohydrodynamic model predicts spikes in the horizontal current density, and associated resistive heating rate per unit volume Q. The spikes appear as increases by orders of magnitude above background values in neutral line regions (NLRs) of active regions (ARs). The largest spikes typically occur a few hours to a few days prior to M or X flares. The spikes correspond to large vertical derivatives of the horizontal magnetic field. The model takes as input the photospheric magnetic field observed by the Helioseismic & Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO) satellite. This 2.5 D field is used to determine an analytic expression for a 3 D magnetic field, from which the current density, vector potential, and electric field are computed in every AR pixel for 14 ARs. The field is not assumed to be force-free. The spurious 6, 12, and 24 hour Doppler periods due to SDO orbital motion are filtered out of the time series of the HMI magnetic field for each pixel using a band pass filter. The subset of spikes analyzed at the pixel level are found to occur on HMI and granulation scales of 1 arcsec and 12 minutes. Spikes are found in ARs with and without M or X flares, and outside as well as inside NLRs, but the largest spikes are localized in the NLRs of ARs with M or X flares. The energy to drive the heating associated with the largest current spikes comes from bulk flow kinetic energy, not the electromagnetic field, and the current density is highly non-force free. The results suggest that, in combination with the model, HMI is revealing strong, convection driven, non-force free heating events on granulation scales, and that it is plausible these events are correlated with subsequent M or X flares. More and longer time series need to be analyzed to determine if such a correlation exists. Above an AR dependent threshold value of Q, the number of events N(Q) with heating rates greater than or equal to Q obeys a scale invariant power law distribution for each AR given by N(Q) varies Q(sup -s), where 0.40 less than or equal to S less than or equal to 0.53, with a mean and standard deviation across the 14 ARs of 0.47 and 0.045, showing there is little variation of S from one AR to another. These properties of N(Q) are in close agreement with those of the distribution N(E) for the total energy E of solar flares, determined from observations to be N(E) = constant x E(sup -alpha). From observations of nanoflares in the 0.7 to 4 MK range, and from observations of flares in hard X-rays, it is found that 0.51 less than or equal to alpha less than or equal to 0.57, and 0.4 less than or equal to alpha less than or equal to 0.6, respectively (Crosby et al. 1993, Sol. Phys., 143, 275; Aschwanden & Parnell 2002, ApJ, 572, 1048). Observations also show that, as is found here for the exponent S, there is little variation of alpha with AR (Wheatland 2000, ApJ, 532, 1209), indicating N(E) and N(Q) are largely independent of individual properties of ARs such as area, total magnetic flux, and distribution of current density (i.e. non-potentiality). Therefore the power law scaling of the photospheric heating rate Q computed here on granulation scales is essentially identical to that found for coronal observations of flare energies on scales 1-2 orders of magnitude larger. This suggests the physical mechanisms that cause Q and coronal flares are closely related. It seems likely that Q is the signature of a magnetic reconnection process in an energy range and volume orders of magnitude smaller than those of flares. In this context, at least the larger spikes in Q might be signatures of UV photospheric or lower chromospheric bombs in which plasma is heated to temperatures approximately 10(exp -5) K (Peter et al. 2014, Science 346, 1255726; Judge 2015, ApJ, 808, 116). In addition, lattice based avalanche simulations of flare energy release predict 0.4 less than or equal to alpha less than or equal to 0.5, while analytic, fractal-diffusive self-organized criticality models predict 0.4 less than or equal to alpha less than or equal to 0.67, in excellent agreement with observations, and the results presented here (Aschwanden & Parnell 2002, ApJ, 572, 1048; Aschwanden 2012, A&A, 539, A2; Aschwanden 2013, in "Self Organized Criticality Systems"; Aschwanden et al. 2016, SSR, 198, 47).

Goodman, Michael↗
NASA NTRS2007

Laser Sounder for Global Measurement of CO2 Concentrations in the Troposphere from Space

Measurements of tropospheric CO2 abundance with global-coverage, a few hundred km spatial and monthly temporal resolution are needed to quantify processes that regulate CO2 storage by the land and oceans. The Orbiting Carbon Observatory (OCO) is the first space mission focused on atmospheric CO2 for measuring total column CO, and O2 by detecting the spectral absorption in reflected sunlight. The OCO mission is an essential step, and will yield important new information about atmospheric CO2 distributions. However there are unavoidable limitations imposed by its measurement approach. These include best accuracy only during daytime at moderate to high sun angles, interference by cloud and aerosol scattering, and limited signal from CO2 variability in the lower tropospheric CO2 column. We have been developing a new laser-based technique for the remote measurement of the tropospheric CO2 concentrations from orbit. Our initial goal is to demonstrate a lidar technique and instrument technology that will permit measurements of the CO2 column abundance in the lower troposphere from aircraft. Our final goal is to develop a space instrument and mission approach for active measurements of the CO2 mixing ratio at the 1-2 ppmv level. Our technique is much less sensitive to cloud and atmospheric scattering conditions and would allow continuous measurements of CO2 mixing ratio in the lower troposphere from orbit over land and ocean surfaces during day and night. Our approach is to use the 1570nm CO2 band and a 3-channel laser absorption spectrometer (i.e. lidar used an altimeter mode), which continuously measures at nadir from a near polar circular orbit. The approach directs the narrow co-aligned laser beams from the instrument's lasers toward nadir, and measures the energy of the laser echoes reflected from land and water surfaces. It uses several tunable fiber laser transmitters which allowing measurement of the extinction from a single selected CO2 absorption line in the 1570 nm band. This band is free from interference from other gases and has temperature insensitive absorption lines. During the measurement the lasers are tuned on- and off- a selected CO2 line near 1572 nm and a selected O2 line near 768 nm in the Oxygen A band at kHz rates. The lasers use tunable diode seed lasers followed by fiber amplifiers, and have spectral widths much narrower than the gas absorption lines. The receiver uses a 1-m diameter telescope and photon counting detectors and measures the background light and energies of the laser echoes from the surface. The extinction and column densities for the CO2 and O2 gases are estimated from the ratio of the on and offline surface echo via the differential optical absorption technique. Our technique rapidly alternates between several on-line wavelengths set to the sides of the selected gas absorption lines. It exploits the atmospheric pressure broadening of the lines to weight the measurement sensitivity to the atmospheric column below 5 km. This maximizes sensitivity to CO2 in the boundary layer, where variations caused by surface sources and sinks are largest. Simultaneous measurements of O2 column will use an identical approach with an O2 line. Thee laser frequencies are tunable and have narrow (MHz) line widths. In combination with sensitive photon counting detectors these enables much higher spectral resolution and precision than is possible with passive spectrometer. 1aser backscatter profiles are also measured, which permits identifying measurements made to cloud tops and through aerosol layers. The measurement approach using lasers in common-nadir-zenith path allows retrieving CO2 column mixing ratios in the lower troposphere irrespective of sun angle. Pulsed laser signals, time gated receiver and a narrow receiver field-of-view are used to isolate the surface laser echo signals and to exclude photons scattered from clouds and aerosols. Nonetheless, the optical absorption change due to a change of a few ppO2 is small, <1 % which makes achieving the needed measurement sensitivities and stabilities quite challenging. Measurement SNRs and stabilities of >600:1 are needed to estimate CO2 mixing ratio at the 1-2 ppm level. We have calculated characteristics of the technique and have demonstrated aspects of the laser, detector and receiver approaches in th e laboratory We have also measured O2 in an absorption cell, and made C02 measurements over a 400 m long (one way) horizontal path using a sensor breadboard. We will describe these and more details of our approach in the paper.

Abshire, James B.↗