Estimation of absorption and backscattering coefficients from in situ radiometric measurements: theory and validation in case II waters

A model that relates the coefficients of absorption (a) and backscattering (b(b)) to diffuse attenuation (K(d)), radiance reflectance (R(L)), and the mean cosine for downward irradiance (mu(d)) is presented. Radiance transfer simulations are used to verify the physical validity of the model for a wide range of water column conditions. Analysis of thee radiance transfer simulations suggest that absorption and backscattering can be estimated with average errors of 1% and 3%, respectively, if the value of mu(d) is known with depth. If the input data set is restricted to variables that can be derived from measurements of upward radiance (L(u)) and downward irradiance (E(d)), it is necessary to use approximate values of mu(d). Examination of three different approximation schemes for mu(d) shows that the average error for estimating a and b(b) increases to approximately 13%. We tested the model by using measurements of L(u) and E(d) collected from case II waters off the west coast of Scotland. The resulting estimates of a and b(b) were compared with independent in situ measurements of these parameters. Average errors for the data set were of the order of 10% for both absorption and backscattering.     

Analytic remote-sensing optical algorithms requiring simple and practical field parameter inputs

Spectral in-water measurements of downward irradiance (E(d)), upward irradiance (E(u)), and nadir radiance (L(u)) are sufficient to calculate the scalar irradiances E(0), E(0d), and E(0u), the average cosines mu, mu(d), and mu(u), the light absorption coefficient a, the backscattering coefficient b(b), and the so-called f factor that relates to R, a, and b(b). The solar elevation of 42 degrees is a special case in which mu(d) is independent of all variables except solar elevation. The algorithms are valid for solar elevations between 12 degrees and 81 degrees for horizontally stratified clear and turbid deep waters.     

Determination of the fluorescence quantum yield by oceanic phytoplankton in their natural habitat

Sun-stimulated chlorophyll a fluorescence has been measured in situ, within the upward and downward light fields, in oceanic waters with chlorophyll concentrations of 0.04-3 mg m(-3). We combined these signals with phytoplankton absorption spectra to derive the fluorescence quantum yield, phi (number of photons emitted by fluorescence/number of absorbed photons). phi was derived separately from hyperspectral (upward and downward) irradiance measurements (with a LI-COR Instruments spectroradiometer) and from nadir radiance near 683 nm (with a Biospherical Instruments profiler). The contribution of inelastic Raman scattering to the signal in the red band was assessed and subtracted. Raman-corrected phi values derived from the two instruments compared well. Vertical phi profiles were strongly structured, with maximal (5-6%) values at depth, whereas phi was approximately =1% in near-surface waters (measurements made approximately at solar noon). These near-surface values are needed for interpretation of remotely sensed fluorescence signals. This optical study shows that the fluorescence yield of algae in their natural environment can be accurately derived in a nonintrusive way with available instrumentation and adequate protocols.     

Ocean optical source estimation with widely spaced irradiance measurements

A method was developed for determining the spatial distribution of a source from downward and upward irradiance measurements at a single wavelength in seawater of known optical properties. The algorithm uses measurements at two depths located an arbitrary distance apart and solves two nonlinear equations for two parameters that fit a globally exponential or linear source shape. Complex spatially dependent source shapes can be estimated from an irradiance profile by piecing together estimates from neighboring measurement pairs. Numerical tests illustrate the sensitivity of the algorithm to depth, measurement spacing, chlorophyll concentration, sensor noise, and uncertainty in the a priori assumed inherent optical properties. The algorithm works well with widely spaced measurements, moderate sensor noise, and uncertainties in the optical properties regardless of whether the assumed and true profiles are the same shape.     

Ocean inherent optical property estimation from irradiances

A method is evaluated for estimating the absorption coefficient a and the backscattering coefficient b(b) from measurements of the upward and downward irradiances E(u)(z) and E(d)(z). With this method, the reflectance ratio R(z) and the downward diffuse attenuation coefficient K(d)(z) obtained from E(u)(z) and E(d)(z) are used to estimate the inherent optical properties R(infinity) and K(infinity) that are the asymptotic values of R(z) and K(d)(z), respectively. For an assumed scattering phase function beta , there are unique correlations between the values of R(infinity) and K(infinity) and those of a and b(b) that can be derived from the radiative transfer equation. Good estimates of a and the Gordon parameter G = b(b)/(a + b(b)) can be obtained from R(infinity) and K(infinity) if the true scattering phase function is not greatly different from the assumed function. The method works best in deep, homogeneous waters, but can be applied to some cases of stratified waters. To improve performance in shallow waters where bottom effects are important, the deep- and shallow-measurement reflectance models also are developed.     

Irradiance of Horizontal Quartz-Halogen Standard Lamps

Spectral irradiance calibrations often require that irradiance standard lamps be oriented differently than the normal calibration orientation used at the National Institute of Standards and Technology and at other standards laboratories. For example, in solar measurements the instruments are generally upward viewing, requiring horizontal working standards for minimization of irradiance calibration uncertainties. To develop a working standard for use in a solar ultraviolet intercomparison, NIST determined the irradiance of quartz-halogen lamps operating in the horizontal position, rather than in the customary vertical position. An experimental technique was developed which relied upon equivalent lamps with independent mounts for each orientation and a spectroradiometer with an integrating sphere whose entrance port could be rotated 90° to view either lamp position. The results presented here are limited to 1000 W quartz-halogen type lamps at ultraviolet wavelengths from 280 nm to 400 nm. Sources of uncertainty arose from the lamps, the spectroradiometer, and the lamp alignment, and increased the uncertainty in the irradiance of horizontal lamps by less than a factor of two from that of vertical NIST standard lamps. The irradiance of horizontal lamps was less than that of vertical lamps by approximately 6 % at long wavelengths (400 nm) to as much as 12 % at the shortest wavelengths (280 nm). Using the Wien radiation law, this corresponds to color temperature differences of 15.7 K and 21.3 K for lamps with clear and frosted envelopes, respectively.     

Semi-empirical correction algorithm for AC-9 measurements in a coccolithophore bloom

Values for the coefficients of absorption (a) and attenuation (c) obtained from AC-9 measurements in coccolithophore blooms do not provide satisfactory inputs for radiance transfer models. We have therefore modified the standard AC-9 scattering correction algorithm by including an extra term, F(lambda, lambda(r)), which allows for possible wavelength dependence in the scattering phase function. We estimated the magnitude of F(lambda, lambda(r)), which is unity in the standard algorithm, by adjusting the absorption and scattering values in Hydrolight radiance transfer calculations until the depth profiles of downward irradiance (E(d)) and upward radiance (L(u)) matched those measured in situ. The modified algorithm was tested with data from a phytoplankton bloom dominated by the coccolithophore Emiliania huxleyi, which occurred in the western English Channel in May 2001. In this paper, we only have sufficient data to adequately constrain the radiance transfer model in one wave band centered on 488 ma. A single value of F(lambda, lambda(r)) = 1.4 was found to produce satisfactory agreement between modeled and observed profiles at four widely spaced stations within the bloom. Measurements of the ratio of backscattering (b(b)) to total scattering (b) showed significant wavelength dependence at these stations.     

Enhanced solar irradiance across the atmosphere-sea ice interface: a quantitative numerical study

A radiative transfer model for the coupled atmosphere-sea ice system is used to study the change in downward irradiance across the air-ice interface. Computations demonstrate that the downward solar irradiance can be significantly enhanced across the air-ice interface. The enhancement is mainly due to light in the ice that is scattered upward and then totally reflected by the air-ice interface. It depends primarily on the change in the index of refraction across this interface and the optical properties of the sea ice, but also on the direct solar and sky illumination of the interface.     

Measured and modeled radiometric quantities in coastal waters: toward a closure

Accurate radiative transfer modeling in the coupled atmosphere-sea system is increasing in importance for the development of advanced remote-sensing applications. Aiming to quantify the uncertainties in the modeling of coastal water radiometric quantities, we performed a closure experiment to intercompare theoretical and experimental data as a function of wavelength lambda and water depth z. Specifically, the study focused on above-water downward irradiance E(d)(lambda, 0+) and in-water spectral profiles of upward nadir radiance L(u)(lambda, z), upward irradiance E(u)(lambda, z), downward irradiance E(d)(lambda, z), the E(u)(lambda, z)/L(u)(lambda, z) ratio (the nadir Q factor), and the E(u)(lambda, z)/E(d)(lambda, z) ratio (the irradiance reflectance). The theoretical data were produced with the finite-element method radiative transfer code ingesting in situ atmospheric and marine inherent optical properties. The experimental data were taken from a comprehensive coastal shallow-water data set collected in the northern Adriatic Sea. Under various measurement conditions, differences between theoretical and experimental data for the above-water E(d)(lambda, 0+) and subsurface E(d)(lambda, 0-) as well as for the in-water profiles of the nadir Q factor were generally less than 15%. In contrast, the in-water profiles of L(u)(lambda, z), E(d)(lambda, z), E(u)(lambda, z) and of the irradiance reflectance exhibited larger differences [to approximately 60% for L(u)(lambda, z) and E(u)(lambda, z), 30% for E(d)(lambda, z), and 50% for the irradiance reflectance]. These differences showed a high sensitivity to experimental uncertainties in a few input quantities used for the simulations: the seawater absorption coefficient; the hydrosol phase function backscattering probability; and, mainly for clear water, the bottom reflectance.     

Asymptotic optical depths in source-free ocean waters

The depth dependence for which the downward diffuse attenuation coefficient, the upward-to-downward plane irradiance ratio, the vertically upward radiance-to-downward plane irradiance ratio, and the mean cosine of the radiance depend negligibly on the surface incident illumination have been examined. The depths at which these coefficients approach to within a specified percent of their asymptotic values depends significantly on the characteristics of the incident illumination and on the inherent optical properties of the water. This information is useful when solving inverse ocean optics problems with a method for which the radiance is assumed to be approximately in the asymptotic regime.     

Comparison of the ocean inherent optical properties obtained from measurements and inverse modeling

A model developed recently by Loisel and Stramski [Appl. Opt. 39, 3001-3011 (2000)] for estimating the spectral absorption a(lambda), scattering b(lambda), and backscattering b(b)(lambda) coefficients in the upper ocean from the irradiance reflectance just beneath the sea surface R(lambda, z = 0(-)) and the diffuse attenuation of downwelling irradiance within the surface layer ?K(d)(lambda)?(1) is compared with measurements. Field data for this comparison were collected in different areas including off-shore and near-shore waters off southern California and around Europe. The a(lambda) and b(b)(lambda) values predicted by the model in the blue-green spectral region show generally good agreement with measurements that covered a broad range of conditions from clear oligotrophic waters to turbid coastal waters affected by river discharge. The agreement is still good if the model estimates of a(lambda) and b(b)(lambda) are based on R(lambda, z = 0(-)) used as the only input to the model available from measurements [as opposed to both R(lambda, z = 0(-)) and ?K(d)(lambda)?(1) being measured]. This particular mode of operation of the model is relevant to ocean-color remote-sensing applications. In contrast to a(lambda) and b(b)(lambda) the comparison between the modeled and the measured b(lambda) shows large discrepancies. These discrepancies are most likely attributable to significant variations in the scattering phase function of suspended particulate matter, which were not included in the development of the model.     

Experimental Investigation of Buoyant Laminar Jet Diffusion Flames in an Inverted Configuration

Buoyant laminar jet diffusion flames are studied experimentally in an inverted configuration, where gaseous fuel-stream jets vertically downward into air. Flame shape, thermal structure, soot and stability behaviors are obtained until the blowoff limit is reached. By comparing with conventional jet flames, which are established when the fuel jets upward, the effects of buoyancy on laminar diffusion flames are analysed. Downward flame yields larger flame height, although the non-dimensional flame height increases linearly with the Reynolds number at nozzle exit, which is similar to upward flame. Possible reasons for the increased flame height include flow deceleration within downward buoyant flames and presence of more combustion products surrounded the jet stream, thus slowing mixing process between fuel and air. The different relative directions of buoyant flows and jet streams also result in different temperature distributions in downward and upward flames, and a stagnant interface produced by the balance between buoyant flow and jet stream is particularly observed downstream of downward flame. Downward flames contain more soot and the soot formation region is wider, which are mainly attributed to the modifications of flow field and soot path. In addition, downward and upward flames stabilize at different axial positions relative to the nozzle exit. Because of increased characteristic flame residence time, downward flames have higher blowoff limits. The downward jet flame provides an alternative configuration to upward jet flame in studying buoyant diffusion flames due to the different manifestations of buoyancy effects.     

Effects of cosine error in irradiance measurements from field ocean color radiometers

The cosine error of in situ seven-channel radiometers designed to measure the in-air downward irradiance for ocean color applications was investigated in the 412-683 nm spectral range with a sample of three instruments. The interchannel variability of cosine errors showed values generally lower than +/-3% below 50 degrees incidence angle with extreme values of approximately 4-20% (absolute) at 50-80 degrees for the channels at 412 and 443 nm. The intrachannel variability, estimated from the standard deviation of the cosine errors of different sensors for each center wavelength, displayed values generally lower than 2% for incidence angles up to 50 degrees and occasionally increasing up to 6% at 80 degrees. Simulations of total downward irradiance measurements, accounting for average angular responses of the investigated radiometers, were made with an accurate radiative transfer code. The estimated errors showed a significant dependence on wavelength, sun zenith, and aerosol optical thickness. For a clear sky maritime atmosphere, these errors displayed values spectrally varying and generally within +/-3%, with extreme values of approximately 4-10% (absolute) at 40-80 degrees sun zenith for the channels at 412 and 443 nm. Schemes for minimizing the cosine errors have also been proposed and discussed.     

Estimation of the inherent optical properties of natural waters from the irradiance attenuation coefficient and reflectance in the presence of Raman scattering

By means of radiative transfer simulations we developed a model for estimating the absorption a, the scattering b, and the backscattering b(b) coefficients in the upper ocean from irradiance reflectance just beneath the sea surface, R(0-), and the average attenuation coefficient for downwelling irradiance, 1, between the surface and the first attenuation depth. The model accounts for Raman scattering by water, and it does not require any assumption about the spectral shapes of a, b, and b(b). The best estimations are obtained for a and b(b) in the blue and green spectral regions, where errors of a few percent to <10% are expected over a broad range of chlorophyll concentration in water. The model is useful for satellite ocean color applications because the model input, R(0-) and 1, can be retrieved from remote sensing and the model output, a and b(b), is the major determinant of remote-sensing reflectance.     

Absolute spectral measurements of direct solar ultraviolet irradiance with a Brewer spectrophotometer

A methodology for the absolute calibration of spectral measurements of direct solar ultraviolet radiation, performed with a Brewer spectrophotometer is presented. The method uses absolute measurements of global and diffuse solar irradiance obtained practically simultaneously at each wavelength with the direct-Sun component. On the basis of this calibration, direct-Sun spectra, measured over a wide range of solar zenith angles at a high altitude site, were used to determine the extraterrestrial solar spectrum by applying the Langley extrapolation method. Finally this spectrum is compared with a solar spectrum derived from the airborne tunable laser absorption spectrometer 3 Space Shuttle mission, showing an agreement of better than +/-3%.     

新型静电驱动双向平动MEMS变形镜(英文)

设计了一种新型静电驱动双向平动的MEMS变形镜,其结构包括3部分:中央下电极,4个静电驱动杠杆结构以及反射镜面(上电极).该变形镜有两种驱动模式:向下驱动模式和向上驱动模式.在向上驱动时,4个杠杆结构实现了位移的放大;向下驱动时,利用了非线性弹性系数法扩展了镜面的平动范围.采用表面硅工艺完成了变形镜的加工.通过白光干涉仪对变形镜的测试表明:在向上驱动模式下,变形镜在驱动电压为31 V时位移达到1.1μm;在向下驱动模式下,变形镜在6 V的驱动电压下位移达到1.1μm.变形镜的总行程为2.2μm,达到了同样工艺下传统变形镜的3倍左右.     

Traveling reference spectroradiometer for routine quality assurance of spectral solar ultraviolet irradiance measurements

A transportable reference spectroradiometer for measuring spectral solar ultraviolet irradiance has been developed and validated. The expanded uncertainty of solar irradiance measurements with this reference spectroradiometer, based on the described methodology, is 8.8% to 4.6%, depending on the wavelength and the solar zenith angle. The accuracy of the spectroradiometer was validated by repeated site visits to two European UV monitoring sites as well as by regular comparisons with the reference spectroradiometer of the European Reference Centre for UV radiation measurements in Ispra, Italy. The spectral solar irradiance measurements of the Quality Assurance of Spectral Ultraviolet Measurements in Europe through the Development of a Transportable Unit (QASUME) spectroradiometer and these three spectroradiometers have agreed to better than 6% during the ten intercomparison campaigns held from 2002 to 2004. If the differences in irradiance scales of as much as 2% are taken into account, the agreement is of the order of 4% over the wavelength range of 300-400 nm.     

Flow and heat transfer on cryogenic flow boiling during tube quenching under upward and downward flow

The gravity effects on quenching of tube by cryogenic fluids for the development of cryogenic fluid management on orbit are studied. In this paper, the effects of the tube diameter, the flow directions, and the mass velocity on the tube quenching using liquid nitrogen are investigated systematically in the terrestrial conditions. The experiments are performed by the mass velocity between 100–600 kg/m²s in downward and upward flow directions by using three difference inner diameters of the transparent heated tube (7, 10, 13.6 mm) for measuring fluid behavior observations and heat transfer measurements simultaneously. The results indicate that the difference between the minimum heat fluxes under downward and upward flow conditions increased as the mass velocity increased. These characteristics of heat transfer were caused by filamentary flow pattern that was found in only downward flow and high mass velocity conditions.     

Charge-coupled device spectrograph for direct solar irradiance and sky radiance measurements

The characterization of a charged-coupled device (CCD) spectrograph developed at the Laboratory of Atmospheric Physics, Thessaloniki is presented. The absolute sensitivity of the instrument for direct irradiance and sky radiance measurements was determined, respectively, with an uncertainty of 4.4% and 6.6% in the UV-B, and 3% and 6% in the UV-A, visible and near-infrared (NIR) wavelength ranges. The overall uncertainty associated with the direct irradiance and the sky radiance measurements is, respectively, of the order of 5% and 7% in the UV-B, increasing to 10% for low signals [e.g., at solar zenith angles (SZAs) larger than 70 degrees ], and 4% and 6% in the UV-A, visible, and NIR. Direct solar spectral irradiance measurements from an independently calibrated spectroradiometer (Bentham DTM 300) were compared with the corresponding CCD measurements. Their agreement in the wavelength range of 310-500nm is within 0.5% +/- 1.1% (for SZA between 20 degrees and 70 degrees ). Aerosol optical depth (AOD) derived by the two instruments using direct Sun spectra and by a collocated Cimel sunphotometer [Aerosol Robotic network (AERONET)] agree to within 0.02 +/- 0.02 in the range of 315-870 nm. Significant correlation coefficients with a maximum of 0.99 in the range of 340-360 nm and a minimum of 0.90 at 870 nm were found between synchronous AOD measurements with the Bentham and the Cimel instruments.     

Operating characteristics of X-bar charts with asymmetric control limits

When the assumptions behind the Shewhart chart are not met, policies other than the traditional 3-sigma limits may enable speedier and more economical detection of process change. If a process has unequal probabilities of downward and upward shifts, downward shifts of differing magnitude from upward shifts, and/or standard deviations after downward shifts different from those after upward shifts, then increases in the in-control average run length and decreases in the out-of-control average run length are simultaneously achievable with control limits placed at unequal distances from the process mean. The magnitude of these improvements is investigated for several types of process. The results are then extended to median charts for processes with output that is not normally distributed.