Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed state estimation”

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.

At least 433 records · Page 24

Determination of Non-Thermal Velocity Distributions from SERTS Linewidth Observations

Non-thermal velocities obtained from the measurement of coronal Extreme Ultraviolet (EUV) linewidths have been consistently observed in solar EUV spectral observations and have been theorized to result fro m many plausible scenarios including wave motions, turbulence, or magnetic reconnection. Constraining these velocities can provide a physical limit for the available energy resulting from unresolved motions in the corona. We statistically determine a series of non-thermal velocity distributions from linewidth measurements of 390 emission lines from a wide array of elements and ionization states observed during the Solar Extreme Ultraviolet Research Telescope and Spectrograph 1991-1997 flights covering the spectral range 174-418 A and a temperature range from 80,000 K to 12.6 MK. This sample includes 248 lines from active regions, 101 lines from quiet-Sun regions, and 41 lines were observed from plasma off the solar limb. We find a strongly peaked distribution corresponding to a non-thermal velocity of 19-22 km/s in all three of the quiet-Sun, active region, and off-limb distributions. For the possibility of Alfven wave resonance heating, we find thai velocities in the core of these distributions do not provide sufficient energy, given typical densities and magnetic field strengths for the coronal plasma, to overcome the estimated coronal energy losses required to maintain the corona at the typical temperatures working as the sole mechanism. We find that at perfect efficiency 50%-60% of the needed energy flux can be produced from the non-thermal velocities measured.

Coyner, Aaron J.↗

Impact of Surface Roughness and Soil Texture on Mineral Dust Emission Fluxes Modeling

Dust production models (DPM) used to estimate vertical fluxes of mineral dust aerosols over arid regions need accurate data on soil and surface properties. The Laboratoire Inter-Universitaire des Systemes Atmospheriques (LISA) data set was developed for Northern Africa, the Middle East, and East Asia. This regional data set was built through dedicated field campaigns and include, among others, the aerodynamic roughness length, the smooth roughness length of the erodible fraction of the surface, and the dry (undisturbed) soil size distribution. Recently, satellite-derived roughness length and high-resolution soil texture data sets at the global scale have emerged and provide the opportunity for the use of advanced schemes in global models. This paper analyzes the behavior of the ERS satellite-derived global roughness length and the State Soil Geographic data base-Food and Agriculture Organization of the United Nations (STATSGO-FAO) soil texture data set (based on wet techniques) using an advanced DPM in comparison to the LISA data set over Northern Africa and the Middle East. We explore the sensitivity of the drag partition scheme (a critical component of the DPM) and of the dust vertical fluxes (intensity and spatial patterns) to the roughness length and soil texture data sets. We also compare the use of the drag partition scheme to a widely used preferential source approach in global models. Idealized experiments with prescribed wind speeds show that the ERS and STATSGO-FAO data sets provide realistic spatial patterns of dust emission and friction velocity thresholds in the region. Finally, we evaluate a dust transport model for the period of March to July 2011 with observed aerosol optical depths from Aerosol Robotic Network sites. Results show that ERS and STATSGO-FAO provide realistic simulations in the region.

textures↗

Integrating Cloud-Based Workflows in Continental-Scale Cropland Extent Classification

Accurate information on cropland spatial distribution is required for global-scale assessments and agricultural land use policies. Cloud computing platforms such as Google Earth Engine (GEE) provide unprecedented opportunities for large-scale classifications of Landsat data. We developed a novel method to fuse pixel-based random forest classification of continental-scale Landsat data on GEE and an object-based segmentation approach known as recursive hierarchical segmentation (RHSeg). Using our fusion method, we produced a continental-scale cropland extent map for North America at 30m spatial resolution for the nominal year 2010. The total cropland area for North America was estimated at 275.18 million hectares (Mha). The overall accuracies of the map are>90% across the continent. This map also compares well with the United States Department of Agriculture (USDA) cropland data layer (CDL), Agriculture and Agri-food Canada (AAFC) annual crop inventory (ACI), and the Mexican government agency Servicio de Informacion Agroalimentaria y Pesquera (SIAP)'s agricultural boundaries. Furthermore, our map compared well with sub-country statistics including state-wise and county-wise cropland statistics in regression models resulting in R2 > 0.84. This key contribution paves the way for more detailed products such as crop intensity, crop type, and crop irrigation, and provides a method for creating high-resolution cropland extent maps for other countries where spatial information about croplands are not as prevalent.

Massey, Richard↗

Dust and gas absorption in the High Mass X-ray Binary IGR J16318−4848

Context. With an absorption column density on the order of 10(exp 24)per sq. cm, IGR J16318-4848 is one of the most extreme cases of a highly obscured High Mass X-ray Binary. Besides the overall continuum absorption, the source spectrum exhibits a strong iron and nickel fluorescence line complex at 6.4 keV. Previous empirical modeling of these features and comparison with radiative transfer simulations raised questions about the structure and covering fraction of the absorber and the profile of the fluorescence lines. Aims. We aim at a self-consistent description of the continuum absorption, the absorption edges and the fluorescent lines to constrain properties of the absorbing material, such as ionization structure and geometry. We further investigate the effects of dust absorption on the observed spectra and the possibility of fluorescent emission from dust grains. Methods. We use XMM-Newton and NuSTAR spectra to first constrain empirically the incident continuum and fluorescence lines. Next, we use XSTAR to construct a customized photoionization model where we vary the ionization parameter, column density, and covering fraction. In a third step, we model the absorption and fluorescence in a dusty olivine absorber and employ both, a simple analytical model for the fluorescent line emission and Monte Carlo radiative transfer, spectral shapes and line fluxes that are very close to the data are generated. Results. Our empirical spectral modeling is in agreement with previous works. Our second model, the single gas absorber does not describe the observational data. In particular, irrespective of the ionization state or column density of the absorber, a much higher covering fraction than previously estimated is needed to produce the strong fluorescence lines and the large continuum absorption. A dusty, spherical absorber (modeled as consisting of olivine dust, although the nature of dust cannot be constrained) is well able to produce the observed continuum absorption and edges. Conclusions. A dense, dusty absorber in the direct vicinity of the source consisting of dust offers a consistent description of both the strong continuum absorption and the strong emission features in the X-ray spectrum of IGR J16318􀀀4848. In particular, for low optical depth of individual grains the dust will contribute significantly to the fluorescent emission, which is the case for typical densities and grain size distribution models.

Ralf Ballhausen↗

Intercomparison of Satellite Derived Gravity Time Series with Inferred Gravity Time Series from TOPEX/POSEIDON Sea Surface Heights and Climatological Model Output

The upcoming GRACE mission promises to open a window on details of the global mass budget that will have remarkable clarity, but it will not directly answer the question of what the state of the Earth's mass budget is over the critical last quarter of the 20th century. To address that problem we must draw upon existing technologies such as SLR, DORIS, and GPS, and climate modeling runs in order to improve our understanding. Analysis of long-period geopotential changes based on SLR and DORIS tracking has shown that addition of post 1996 satellite tracking data has a significant impact on the recovered zonal rates and long-period tides. Interannual effects such as those causing the post 1996 anomalies must be better characterized before refined estimates of the decadal period changes in the geopotential can be derived from the historical database of satellite tracking. A possible cause of this anomaly is variations in ocean mass distribution, perhaps associated with the recent large El Nino/La Nina. In this study, a low-degree spherical harmonic gravity time series derived from satellite tracking is compared with a TOPEX/POSEIDON-derived sea surface height time series. Corrections for atmospheric mass effects, continental hydrology, snowfall accumulation, and ocean steric model predictions will be considered.

Cox, C.↗

Self-Nulling Lock-in Detection Electronics for Capacitance Probe Electrometer

A multi-channel electrometer voltmeter that employs self-nulling lock-in detection electronics in conjunction with a mechanical resonator with noncontact voltage sensing electrodes has been developed for space-based measurement of an Internal Electrostatic Discharge Monitor (IESDM). The IESDM is new sensor technology targeted for integration into a Space Environmental Monitor (SEM) subsystem used for the characterization and monitoring of deep dielectric charging on spacecraft. Use of an AC-coupled lock-in amplifier with closed-loop sense-signal nulling via generation of an active guard-driving feedback voltage provides the resolution, accuracy, linearity and stability needed for long-term space-based measurement of the IESDM. This implementation relies on adjusting the feedback voltage to drive the sense current received from the resonator s variable-capacitance-probe voltage transducer to approximately zero, as limited by the signal-to-noise performance of the loop electronics. The magnitude of the sense current is proportional to the difference between the input voltage being measured and the feedback voltage, which matches the input voltage when the sense current is zero. High signal-to-noise-ratio (SNR) is achieved by synchronous detection of the sense signal using the correlated reference signal derived from the oscillator circuit that drives the mechanical resonator. The magnitude of the feedback voltage, while the loop is in a settled state with essentially zero sense current, is an accurate estimate of the input voltage being measured. This technique has many beneficial attributes including immunity to drift, high linearity, high SNR from synchronous detection of a single-frequency carrier selected to avoid potentially noisy 1/f low-frequency spectrum of the signal-chain electronics, and high accuracy provided through the benefits of a driven shield encasing the capacitance- probe transducer and guarded input triaxial lead-in. Measurements obtained from a 2- channel prototype electrometer have demonstrated good accuracy (|error| < 0.2 V) and high stability. Twenty-four-hour tests have been performed with virtually no drift. Additionally, 5,500 repeated one-second measurements of 100 V input were shown to be approximately normally distributed with a standard deviation of 140 mV.

Blaes, Brent R.↗

Aerodynamic Stability and Performance of Next-Generation Parachutes for Mars Descent

The Low Density Supersonic Decelerator Project is developing a next-generation supersonic parachute for use on future Mars missions. In order to determine the new parachute configuration, a wind tunnel test was conducted at the National Full-scale Aerodynamics Complex 80- by 120-foot Wind Tunnel at the NASA Ames Research Center. The goal of the wind tunnel test was to quantitatively determine the aerodynamic stability and performance of various canopy configurations in order to help select the design to be flown on the Supersonic Flight Dynamics tests. Parachute configurations included the diskgap- band, ringsail, and ringsail-variant designs referred to as a disksail and starsail. During the wind tunnel test, digital cameras captured synchronized image streams of the parachute from three directions. Stereo hotogrammetric processing was performed on the image data to track the position of the vent of the canopy throughout each run. The position data were processed to determine the geometric angular history of the parachute, which were then used to calculate the total angle of attack and its derivatives at each instant in time. Static and dynamic moment coefficients were extracted from these data using a parameter estimation method involving the one-dimensional equation of motion for a rotation of parachute. The coefficients were calculated over all of the available canopy states to reconstruct moment coefficient curves as a function of total angle of attack. From the stability curves, useful metrics such as the trim total angle of attack and pitch stiffness at the trim angle could be determined. These stability metrics were assessed in the context of the parachute's drag load and geometric porosity. While there was generally an inverse relationship between the drag load and the stability of the canopy, the data showed that it was possible to obtain similar stability properties as the disk-gap-band with slightly higher drag loads by appropriately tailoring the geometric porosity distribution.

Low Density Supersonic Decelerator (LDSD)↗

A Joint Analysis of GPS Displacement and GRACE Geopotential Data for Simultaneous Estimation of Geocenter Motion and Gravitational Field

Gravitational potential data from GRACE are being used to study mass redistribution within and between the atmosphere, hydrosphere, cryosphere, and solid Earth. The GRACE data are made available in a reference frame with its origin at the center of mass of the Earth system (geocenter) while many other geophysical models and data sets refer to a reference frame attached to the Earth's surface. Changes in the offset between these reference frames (geocenter motion) must be accounted for when GRACE data are used to quantify surface mass changes. In this study, we developed a technique for co‐estimation of geocenter motion and gravitational potential field seamlessly from degree 1 to 90 by simultaneously inverting a set of globally‐distributed GPS displacement time series and the temporally‐varying GRACE gravity data. We found that the effect of geocenter motion was evident particularly in the GPS time series of horizontal displacements. Our estimates of geocenter motion are most consistent with the Satellite Laser Ranging (SLR) results within 1 mm in X and Z components and a submillimeter in Y component, when compared to monthly variability averaged over the period of 2003–2016. The overall magnitude of the degree‐1 (l = 1) surface mass load is estimated to be ~3 cm in equivalent water height annually migrating south‐westward from Europe (December–January) to the South Pacific (June–July). Our results also show that dense GPS network data improve water storage recovery in major river basins in the United States and Europe by contributing significantly to the recovery of higher‐degree (l ≥ ~20) geopotential coefficients.

Razeghi, Mahdiyeh↗

PPDIST, global 0.1° daily and 3-hourly precipitation probability distribution climatologies for 1979–2018

We introduce the Precipitation Probability DISTribution (PPDIST) dataset, a collection of global high-resolution (0.1°) observation-based climatologies (1979–2018) of the occurrence and peak intensity of precipitation (P) at daily and 3-hourly time-scales. The climatologies were produced using neural networks trained with daily P observations from 93,138 gauges and hourly P observations (resampled to 3-hourly) from 11,881 gauges worldwide. Mean validation coefficient of determination (R^(2)) values ranged from 0.76 to 0.80 for the daily P occurrence indices, and from 0.44 to 0.84 for the daily peak P intensity indices. The neural networks performed significantly better than current state-of-the-art reanalysis (ERA5) and satellite (IMERG) products for all P indices. Using a 0.1 mm 3 per h threshold, P was estimated to occur 12.2%, 7.4%, and 14.3% of the time, on average, over the global, land, and ocean domains, respectively. The highest P intensities were found over parts of Central America, India, and Southeast Asia, along the western equatorial coast of Africa, and in the intertropical convergence zone.

Hylke E. Beck↗

An Intercomparison of Microphysical Retrieval Algorithms for Upper Tropospheric Ice Clouds

The large horizontal extent, location in the cold upper troposphere, and ice composition make cirrus clouds important modulators of the earth's radiation budget and climate. Cirrus cloud microphysical properties are difficult to measure and model because they are inhomogeneous in nature and their ice crystal size distribution and habit are not well characterized. Accurate retrievals of cloud properties are crucial for improving the representation of cloud scale processes in large-scale models and for accurately predicting the earth's future climate. A number of passive and active remote sensing retrieval algorithms exist for estimating the microphysical properties of upper tropospheric clouds. We believe significant progress has been made in the evolution of these retrieval algorithms in the last decade, however, there is room for improvement. Members of the Atmospheric Radiation measurement program (ARM) Cloud properties Working Group are involved in an intercomparison of optical depth(tau), ice water path, and characteristic particle size in clouds retrieved using ground-based instruments. The goals of this intercomparison are to evaluate the accuracy of state-of-the-art algorithms, quantify the uncertainties, and make recommendations for improvement.

Comstock, Jennifer M.↗

Fuel cell environmental life cycle assessment : a progress report

Energy generation using fuel cells offers the possibility of high conversion efficiency and reduced pollution. Thus, in the futue, conventional methods of energy generation using coalfired power plants, natural gas turbines and internal combustion engines can potentially be replaced by fuel cells. The United States, Japan and several countries in Europe are currently developing fuel cells for transportation, distributed energy generation and portable electronics applications. On a system basis, the efficiency of chemical to electrical energy conversion in fuel cells can be as high as 45%. Also, because of the lower operating temperatures compared to turbines and internal combustion engines, emissions such as carbon monoxide and nitrogen oxides are extremely small. An assessment of the environmental impact of fuel cells over the entire lifecycle is needed to determine the specific benefits in various application scenarios. Such a comparative pollutants assessment must include a "cradle to grave" approach encompassing fuel production, fuel transportation, materials/manufacturing processes, waste management, and recycling. A comparative pollutant model would be a useful tool for industry and government in estimating energy savings and environmental credits, and providing a holistic understanding on regulatory needs. Another important benefit would be the ability to make appropriate comparisons with competitive energy generation technologies such as hybrids, turbines, and internal combustion engines, and to set up realistic targets. Also, exercising the model will assist in the initial screening for potential hot spots and showstoppers and focus on the right component technologies design for manufacturing, designs for input substitution and waste minimization, and prevent any unforeseen burdens.

Narayanan, S. R.↗

Wind Variability of B Supergiants: The Rapid Rotator HD 64760 (B0.5 Ib) - No. 1

We present the results of a 6 day time series of observations of the rapidly rotating B0.5 Ib star HD 64760. We point out several reasons why such intermediate luminosity B supergiants are ideal targets for wind variability studies and then present our results that show the following: continuous wind activity throughout the 6 day run with the wind never in steady state for more than a few hr; wind variability very near nu = 0 km sec(exp -1) in the resonance lines from the lower ionization stages (Al III and C II); a distinct correlation between variability in the Si III ; lambda(lambda)1300 triplets, the strong C III (lambda)1247 singlet, and the onset of extremely strong wind activity, suggesting a connection between photospheric and wind activity; long temporal coherence in the behavior of the strong absorption events; evidence for large-scale spatial coherence, implied by a whole scale, simultaneous weakening in the wind absorption over a wide range in velocities; and ionization variability in the wind accompanying the largest changes in the absorption strengths of the wind lines. In addition, modeling of the wind lines provides the following information about the state the wind in HD 64760. The number of structures on the portion of a constant velocity surface occulting the stellar disk at a particular time must be quite small, while the number on the entire constant velocity surface throughout the wind must be large. The escape probability at low velocity is overestimated by a normal beta approx. 1 velocity law, perhaps due to the presence of low-velocity shocks deep in the wind or a shallow velocity gradient at low velocity. Estimates of the ionization structure in the wind indicate that the ionization ratios are not those expected from thermal equilibrium wind models or from an extrapolation of previous O star results. The large observed q(N V)/q(Si IV) ratio is almost certainly due to distributed X-rays, but the level of ionization predicted by distributed X-ray wind models is inconsistent with the predicted mass-loss rate. Thus, it is impossible to reconcile the observed ionization ratios and the predicted mass-loss rate within the framework of the available models.

Massa, Derck↗

An Experimental Investigation Into the Temperature Profile of a Compliant Foil Air Bearing

A series of tests was performed to determine the internal temperature profile in a compliant bump-type foil journal air bearing operating at room temperature under various speeds and load conditions. The temperature profile was collected by instrumenting a foil bearing with nine, type K thermocouples arranged in the center and along the bearing s edges in order to measure local temperatures and estimate thermal gradients in the axial and circumferential directions. To facilitate the measurement of maximum temperatures from viscous shearing in the air film, the thermocouples were tack welded to the backside of the bumps that were in direct contact with the top foil. The mating journal was coated with a high temperature solid lubricant that, together with the bearing, underwent high temperature start-stop cycles to produce a smooth, steady-state run-in surface. Tests were conducted at speeds from 20 to 50 krpm and loads ranging from 9 to 222 N. The results indicate that, over the conditions tested, both journal rotational speed and radial load are responsible for heat generation with speed playing a more significant role in the magnitude of the temperatures. The temperature distribution was nearly symmetric about the bearing center at 20 and 30 krpm but became slightly skewed toward one side at 40 and 50 krpm. Surprisingly, the maximum temperatures did not occur at the bearing edge where the minimum film thickness is expected but rather in the middle of the bearing where analytical investigations have predicted the air film to be much thicker. Thermal gradients were common during testing and were strongest in the axial direction from the middle of the bearing to its edges, reaching 3.78 8C/mm. The temperature profile indicated the circumferential thermal gradients were negligible.

Radil, Kevin↗

SO2 Emissions and Lifetimes: Estimates from Inverse Modeling Using In Situ and Global, Space-Based (SCIAMACHY and OMI) Observations

Top-down constraints on global sulfur dioxide (SO2) emissions are inferred through inverse modeling using SO2 column observations from two satellite instruments (SCIAMACHY and OMI). We first evaluated the S02 column observations with surface SO2 measurements by applying local scaling factors from a global chemical transport model (GEOS-Chem) to SO2 columns retrieved from the satellite instruments. The resulting annual mean surface SO2 mixing ratios for 2006 exhibit a significant spatial correlation (r=0.86, slope=0.91 for SCIAMACHY and r=0.80, slope = 0.79 for OMI) with coincident in situ measurements from monitoring networks throughout the United States and Canada. We evaluate the GEOS-Chem simulation of the SO2 lifetime with that inferred from in situ measurements to verity the applicability of GEOS-Chem for inversion of SO2 columns to emissions. The seasonal mean SO2 lifetime calculated with the GEOS-Chem model over the eastern United States is 13 h in summer and 48 h in winter, compared to lifetimes inferred from in situ measurements of 19 +/- 7 h in summer and 58 +/- 20 h in winter. We apply SO2 columns from SCIAMACHY and OMI to derive a top-down anthropogenic SO2 emission inventory over land by using the local GEOS-Chem relationship between SO2 columns and emissions. There is little seasonal variation in the top-down emissions (<15%) over most major industrial regions providing some confidence in the method. Our global estimate for annual land surface anthropogenic SO2 emissions (52.4 Tg S/yr from SCIAMACHY and 49.9 Tg S / yr from OMI) closely agrees with the bottom-up emissions (54.6 Tg S/yr) in the GEOS-Chem model and exhibits consistency in global distributions with the bottom-up emissions (r = 0.78 for SCIAMACHY, and r = 0.77 for OMI). However, there are significant regional differences.

Lee, Chulkyu↗

Mapping the Extent of Mangrove Ecosystem Degradation by Integrating an Ecological Conceptual Model with Satellite Data

Anthropogenic and natural disturbances can cause degradation of ecosystems, reducing their capacity to sustain biodiversity and provide ecosystem services. Understanding the extent of ecosystem degradation is critical for estimating risks to ecosystems, yet there are few existing methods to map degradation at the ecosystem scale and none using freely available satellite data for mangrove ecosystems. In this study, we developed a quantitative classification model of mangrove ecosystem degradation using freely available earth observation data. Crucially, a conceptual model of mangrove ecosystem degradation was established to identify suitable remote sensing variables that support the quantitative classification model, bridging the gap between satellite-derived variables and ecosystem degradation with explicit ecological links. We applied our degradation model to two case-studies, the mangroves of Rakhine State, Myanmar, which are severely threatened by anthropogenic disturbances, and Shark River within the Everglades National Park, USA, which is periodically disturbed by severe tropical storms. Our model suggested that 40% (597 km2) of the extent of mangroves in Rakhine showed evidence of degradation. In the Everglades, the model suggested that the extent of degraded mangrove forest increased from 5.1% to 97.4% following the Category 4 Hurricane Irma in 2017. Quantitative accuracy assessments indicated the model achieved overall accuracies of 77.6% and 79.1% for the Rakhine and the Everglades, respectively. We highlight that using an ecological conceptual model as the basis for building quantitative classification models to estimate the extent of ecosystem degradation ensures the ecological relevance of the classification models. Our developed method enables researchers to move beyond only mapping ecosystem distribution to condition and degradation as well. These results can help support ecosystem risk assessments, natural capital accounting, and restoration planning and provide quantitative estimates of ecosystem degradation for new global biodiversity targets.

Calvin K. F. Lee↗

Completing the Feedback Loop: The Impact of Chlorophyll Data Assimilation on the Ocean State

In anticipation of the integration of a full biochemical model into the next generation GMAO coupled system, an intermediate solution has been implemented to estimate the penetration depth (1Kd_PAR) of ocean radiation based on the chlorophyll concentration. The chlorophyll is modeled as a tracer with sources-sinks coming from the assimilation of MODIS chlorophyll data. Two experiments were conducted with the coupled ocean-atmosphere model. In the first, climatological values of Kpar were used. In the second, retrieved daily chlorophyll concentrations were assimilated and Kd_PAR was derived according to Morel et al (2007). No other data was assimilated to isolate the effects of the time-evolving chlorophyll field. The daily MODIS Kd_PAR product was used to validate the skill of the penetration depth estimation and the MERRA-OCEAN re-analysis was used as a benchmark to study the sensitivity of the upper ocean heat content and vertical temperature distribution to the chlorophyll input. In the experiment with daily chlorophyll data assimilation, the penetration depth was estimated more accurately, especially in the tropics. As a result, the temperature bias of the model was reduced. A notably robust albeit small (2-5 percent) improvement was found across the equatorial Pacific ocean, which is a critical region for seasonal to inter-annual prediction.

MERRA-OCEAN↗

Seven-Year SSM/I-Derived Global Ocean Surface Turbulent Fluxes

A 7.5-year (July 1987-December 1994) dataset of daily surface specific humidity and turbulent fluxes (momentum, latent heat, and sensible heat) over global oceans has been retrieved from the Special Sensor Microwave/Imager (SSM/I) data and other data. It has a spatial resolution of 2.0 deg.x 2.5 deg. latitude-longitude. The retrieved surface specific humidity is generally accurate over global oceans as validated against the collocated radiosonde observations. The retrieved daily wind stresses and latent heat fluxes show useful accuracy as verified by those measured by the RV Moana Wave and IMET buoy in the western equatorial Pacific. The derived turbulent fluxes and input variables are also found to agree generally with the global distributions of annual-and seasonal-means of those based on 4-year (1990-93) comprehensive ocean-atmosphere data set (COADS) with adjustment in wind speeds and other climatological studies. The COADS has collected the most complete surface marine observations, mainly from merchant ships. However, ship measurements generally have poor accuracy, and variable spatial coverages. Significant differences between the retrieved and COADS-based are found in some areas of the tropical and southern extratropical oceans, reflecting the paucity of ship observations outside the northern extratropical oceans. Averaged over the global oceans, the retrieved wind stress is smaller but the latent heat flux is larger than those based on COADS. The former is suggested to be mainly due to overestimation of the adjusted ship-estimated wind speeds (depending on sea states), while the latter is suggested to be mainly due to overestimation of ship-measured dew point temperatures. The study suggests that the SSM/I-derived turbulent fluxes can be used for climate studies and coupled model validations.

Chou, Shu-Hsien↗

The Ground Flash Fraction Retrieval Algorithm Employing Differential Evolution: Simulations and Applications

The ability to estimate the fraction of ground flashes in a set of flashes observed by a satellite lightning imager, such as the future GOES-R Geostationary Lightning Mapper (GLM), would likely improve operational and scientific applications (e.g., severe weather warnings, lightning nitrogen oxides studies, and global electric circuit analyses). A Bayesian inversion method, called the Ground Flash Fraction Retrieval Algorithm (GoFFRA), was recently developed for estimating the ground flash fraction. The method uses a constrained mixed exponential distribution model to describe a particular lightning optical measurement called the Maximum Group Area (MGA). To obtain the optimum model parameters (one of which is the desired ground flash fraction), a scalar function must be minimized. This minimization is difficult because of two problems: (1) Label Switching (LS), and (2) Parameter Identity Theft (PIT). The LS problem is well known in the literature on mixed exponential distributions, and the PIT problem was discovered in this study. Each problem occurs when one allows the numerical minimizer to freely roam through the parameter search space; this allows certain solution parameters to interchange roles which leads to fundamental ambiguities, and solution error. A major accomplishment of this study is that we have employed a state-of-the-art genetic-based global optimization algorithm called Differential Evolution (DE) that constrains the parameter search in such a way as to remove both the LS and PIT problems. To test the performance of the GoFFRA when DE is employed, we applied it to analyze simulated MGA datasets that we generated from known mixed exponential distributions. Moreover, we evaluated the GoFFRA/DE method by applying it to analyze actual MGAs derived from low-Earth orbiting lightning imaging sensor data; the actual MGA data were classified as either ground or cloud flash MGAs using National Lightning Detection Network[TM] (NLDN) data. Solution error plots are provided for both the simulations and actual data analyses.

Koshak, William↗