Oxygen isotope measurements of phosphate from fish teeth and bones

In situ measurements of lunar surface brightness temperatures made as a part of the Apollo Lunar Surface Experiments Package at the Apollo 15 Hadley Rille landing site are reported. Data derived from 5 thermocouples of the Heat Flow Experiment, which are lying on or just above the surface, are used to examine the thermal properties of the upper 15 cm of the lunar regolith using eclipse and nighttime cool-down temperatures. Application of finite-difference techniques in modeling the lunar soil shows the thermocouple data are best fit by a model consisting of a low-density and low-thermal conductivity surface layer approximately 2 cm thick overlying a region increasing in conductivity and density with depth. Conductivities on the order of 1 × 10^?5 W/cm-°K are postulated for the upper layer, with conductivity increasing to the order of 1 × 10^?4 W/cm-°K at depths exceeding 20 cm. An increase in mean temperature with depth indicates that the ratio of radiative to conductive transfer at 350°K is 2.7 for at least the upper few centimeters of lunar soil; this value is nearly twice that measured for returned lunar fines. The thermal properties model deduced from Apollo 15 surface temperatures is consistent with earth-based microwave observations if electrical properties measured on returned lunar fines are assumed.     

月球表层及月壳物质密度分布特征

月球表层与月壳岩石密度的横向与径向的变化,反映了月表及内部成分以及月球演化等特征.本文利用月球勘探者号伽马射线谱仪探测的月表Fe, Th与Mg元素分布数据,依据前人给出的元素含量与岩石类型的关系,对月球表层进行了岩性填图,并结合岩石样品与陨石的密度测试数据建立初始密度模型,采用铁元素与岩石密度的关系对其进行修正,从而建立了月表物质密度分布模型.基于嫦娥一号激光测高数据和日本SELENE计划发布的月球重力模型,计算出月球布格重力异常,进而反演得到月壳0~40 km深度范围内岩石平均密度分布模型.分析表明,大部分区域上,月壳至少月壳上部岩石成分主要以轻质的富含铝、钙、镁质的硅酸盐类岩石为主.由此推测,原始月壳极有可能是由轻质的、富含钙、镁质硅酸盐类岩石构成的全球性月壳.现今的玄武岩与克里普岩只是覆盖于原始的月壳之上的岩层,且厚度不大.     

A Monte Carlo model for the exposure history of lunar dust grains in the ancient solar wind

The theoretical motion of individual dust grains in the lunar regolith is analyzed by using a Monte Carlo statistical code where the variables are the mass and speed distribution of meteorites at the lunar surface and the geometrical shape of impact craters. From these computations the detailed irradiation history of the grains in the ancient solar wind is traced back, over a period of 4 billion years, as a function of the grain size. Then by combining this irradiation scheme with the result of solar wind simulation experiments, the time and depth dependent accumulation of solar wind effects in the theoretical grains (solar wind maturation) is inferred. Finally, the validity of these predictions is tentatively checked by discussing a variety of physical and chemical solar wind effects which are registered in the surface layers of lunar dust grains. Therefore these studies give a tentative scenario for the “maturation” of the lunar regolith with respect to solar wind effects, but they also reveal useful guidelines to deduce meaningful information from such effects. In particular, they suggest a “lunar skin” sampling technique for extracting dust grains in lunar core tubes which could help in deciphering the past activity of the ancient solar wind over a time scale of several billion years.     

Influence of (FeO+TiO2) abundance on the microwave thermal emissions of lunar regolith

One of the essential controls on the microwave thermal emissions (MTE) of the lunar regolith is the abundance of FeO and TiO₂, known as the (FeO+TiO₂) abundance (FTA). In this paper, a radiative transfer simulation is employed first to study the change in the brightness temperature (T_B) with FTA under a range of frequencies and surface temperatures. Then, we analyze the influence of FTA on the MTE of the lunar regolith using microwave sounder (CELMS) data from the Chang’E-2 lunar orbiter, Clementine UV-VIS data, and lunar samples recovered from the Apollo and Surveyor projects. We conclude that: (1) FTA strongly influences the MTE of the lunar regolith, but it is not the decisive control, and (2) FTA decreases slightly with depth. This research plays an essential role in appropriately inverting CELMS data to obtain lunar regolith parameters.     

中国大陆地表温度年变基准场研究

非构造活动或非地震因素对地表热辐射场（地表温度）的影响,对于利用卫星热红外遥感探索地震前兆抑或断层活动有着重要的现实意义.地表温度中,典型非构造活动或非地震因素成份有：由太阳辐射引起的稳定年周期成份和与地形、纬度及能量平衡等因素有关的长期稳定成份,合称年变基准场.本文根据2000～2008年的MODIS/Terra地表温度产品,利用小波分析提取了中国大陆地表温度的年变基准场.在此基础上,结合热传导方程和数学物理方法,获得了年变基准场的（半）定量化表达式.进一步, 利用地表温度的长期稳定成份,获得了温度与海拔、纬度变化之间的定量关系：海拔每增加100 m,温度降低0.51±0.01 K; 纬度每增加1°,温度降低0.77±0.08 K.总之,年变基准场可为利用热红外辐射提取地壳活动信息提供一种参考背景, 温度与海拔、纬度变化之间的定量关系则可用来校正地形起伏和纬度变化对地表温度的影响.     

Influence of the topography of the surface of the Earth on vertical sounding of the temperature profile

Atmospheric temperature and humidity fields as well as information on other meteorological parameters are nowadays retrieved from radiance measurements recorded by operational meteorological satellites. Up to now, the inversion procedures used only take into account crude information on the topography of the Earth’s surface. However, the applied radiative transfer codes have to consider the Earth’s surface as the lower boundary of the atmospheric model and, therefore, need a more precise mean elevation and a classification of the roughness of the Earth’s surface. The influence of the topography of the Earth surface on retrieved temperature profiles is studied by using a physico-statistical inversion method. An objective analysis is made of the more precise mean elevation and derivation of roughness parameters using a new high-resolution digital elevation model (DEM) with a resolution of 500 m×500 m. By means of a geomorphological process and a newly developed topography rejection test, areas with a high surface roughness are localized and singled out. The influence of topography on the retrieved temperature profiles is illustrated by case studies. Changes are found predominantly in areas with a high variation of topography. Using the new high-resolution DEM and the topography rejection test, the geographical position of the calculated temperature profiles tends to be shifted towards areas with a small vertical variation of topography. The mean elevation determined by the new elevation model better characterizes the area observed. Hence, the temperature profiles can be calculated down to lower atmospheric levels. Furthermore, a guess profile better describing the atmospheric situation is selected by the more precise elevation. In addition, the temperature profiles obtained near the coast are improved considerably by the more precise determination of the surface property ‘sea’ and ‘land,’ respectively. Integration of an independent physical information such as topography leads, on average, to a slight improvement of the results of the physico-statistical inversion procedure. In some cases, however, significant improvements have been achieved regarding the desired accuracy of temperature profiles of the order of 1 K. In future, the spatial resolution of new high-resolution sounding instrumentation on the next generation of operational meteorological satellites will be increased. To exploit the resolving power of this new instrumentation, the different variation of the topography of the Earth surface, especially in regions with a high variation of topography, can be taken into account more precisely by using a high-resolution DEM.     

Calibration of an integrated land surface process and radiobrightness (LSP/R) model during summertime

In this study, a soil vegetation and atmosphere transfer (SVAT) model was linked with a microwave emission model to simulate microwave signatures for different terrain during summertime, when the energy and moisture fluxes at the land surface are strong. The integrated model, land surface process/radiobrightness (LSP/R), was forced with weather and initial conditions observed during a field experiment. It simulated the fluxes and brightness temperatures for bare soil and brome grass in the Northern Great Plains. The model estimates of soil temperature and moisture profiles and terrain brightness temperatures were compared with the observed values. Overall, the LSP model provides realistic estimates of soil moisture and temperature profiles to be used with a microwave model. The maximum mean differences and standard deviations between the modeled and the observed temperatures (canopy and soil) were 2.6 K and 6.8 K, respectively; those for the volumetric soil moisture were 0.9% and 1.5%, respectively. Brightness temperatures at 19 GHz matched well with the observations for bare soil, when a rough surface model was incorporated indicating reduced dielectric sensitivity to soil moisture by surface roughness. The brightness temperatures of the brome grass matched well with the observations indicating that a simple emission model was sufficient to simulate accurate brightness temperatures for grass typical of that region and surface roughness was not a significant issue for grass-covered soil at 19 GHz. Such integrated SVAT-microwave models allow for direct assimilation of microwave observations and can also be used to understand sensitivity of microwave signatures to changes in weather forcings and soil conditions for different terrain types.     

Solar-wind induction and lunar conductivity

A summary review of electromagnetic induction driven in the Moon by the interplanetary magnetic field is given. The point of view developed centers on inversion of Fourier transforms of the magnetic field in the free-stream solar wind (forcing function) and the response on the lunar surface measured by Lunar Surface Magnetometers. Conductivity profiles are shown to depend upon the central angle between the magnetometer given by local time and the incident wave-normal direction. The induction excites at least magnetic dipole and quadrupole “radiation”, but any scattered field is confined to the Moon's interior, save for propagation down the cavity where a TE-mode surface wave is generated. Confinement of the induced field on the sunward hemisphere and near the subsolar point is nearly complete, decreasing to the limbs, while in the diamagnetic cavity downstream of the Moon, partial confinement takes place. Both time and spatial multipoles of the induced field are present in the lunar interior complicating inversion into conductivity profiles. Profiles are reviewed and resolution limits are given and compared to those obtained from transient analysis. Finally a qualitative comparison to conductivity in the Earth is given.     

In-situ measurement of the rate of235U fission induced by lunar neutrons

The depth profile of the neutron-induced fission rate of²³⁵U was directly measured to a depth of 350 g/cm² by the Apollo 17 Lunar Neutron Probe Experiment. The fission rate rises sharply from the surface to a broad maximum from 110 to 160 g/cm² and drops off at greater depths. The shape of theoretical depth profile of Lingenfelter et al. fits the measured capture rates well at all depths. The absolute magnitude of the experimental fission rates are (11±17)% lower than those calculated theoretically. The excellent agreement between theory and experiment implies that conclusions drawn previously by interpreting lunar sample data with the theoretical capture rates will not require revision. In particular lunar surface processes, rather than uncertainties in the capture rates, are required to explain the relatively low neutron fluences observed for surface soil samples compared to the fluences expected for a uniformly mixed regolith.     

卫星热红外亮温、气温及地温观测的年变变化对比研究

过对NOAA卫星热红外亮温与野外安装气象观测站接收的气温、地下不同深度地温(0.2, 0.5, 1.5,2 m)进行不同方式的对比研究,分析了卫星热红外亮温、气温、地下不同深度地温的变化特征,探讨了亮温与气温及不同深度地温之间的关系．结果表明:① 卫星热红外亮温观测,由于受天气、云层短周期因素变化影响,曲线呈现高频突跳特征,但按最大值拟合出的亮温曲线有较好的年变变化规律;② 浅层地温受气温及太阳辐射的影响较大,能够体现出日变化,表现出很好的季节变化规律;③ 深层地温年变平稳, 年变变化与季节相关．但与气温相比,表现出滞后效应,且深度越深,滞后时间越长;④ 亮温、气温及深度0.2 m地温三者之间呈现很好的相关性．亮温、气温、0.2 m地温的极值几乎同时段出现,都符合季节变化．分析表明,亮温能够真实地反映地表温度的变化情况,能够为利用卫星热红外亮温提取地震异常信息提供可靠准确证据．本研究结果为理解不同观测属性及其相互关系,以及更好地为地震监测应用提供了基础．      

A method of computing the vertical profile of ozone concentration based on radiation field analysis

Summary A methodical approach is presented of solving the inverse problem of atmospheric optics for the vertical profile of ozone concentration. Observations of spectral sky radiance and direct solar radiation are taken as input data. A gradient method is suggested for solving the inverse problem. on leave from the Astronomical Institute, Slovak Academy of Sciences     

11-Year cyclicity of the near-surface heat flux at the Vostok station

It is known that the solar component is difficult to find in a climatic signal due to its small size and significant internal random disturbances of the climatic system, such as variations in cloudiness, precipitation, winds, and oceanic currents. In the center of the polar cap, these disturbances are as minimal as possible; therefore, solving the inverse heat conduction problem, we can calculate the informative heat flux through the ice cover based on temperature data. Variations in the heat flux are quite significant in amplitude (0.2–0.3 W/m²) and coincide with the 11-year solar cycle in phase. Application of this approach to the indices of the global and sea surface temperatures cannot yield a solar signal, yet it can be traced for solar cycles with big amplitude (cycles 19 and 21). Thus, the found variation in the heat flux is most probably caused by a change in the solar constant (TSI); however, other mechanisms are also possible.     

Synchronous retrieval of land surface temperature and emissivity

This is an old topic for more than ten years to retrieve land surface temperature (LST) from satellite data, but it has not been solved yet. At first, people tried to transplant traditional split window method of sea surface temperature (SST) to the retrieval of LST, but it was found that the emissivities of land surface (ε_i) must be involved in atmospheric correction. Then many different formulas appeared with assumption of emissivities known. In fact, emissivities of land surface with pixel size cannot be known beforehand because of various reasons, so in recent years the focus of attention has been transferred to retrieving emissivities (ε_i) and LST at the same time. Therefore, we have to solve missing equations problem. For this some people try to introduce middle infrared information, but new problems will be brought in which means that it is very difficult to describe middle infrared BRDF of targets with high accuracy and the scattering of atmospheric aerosol cannot be ignored. Therefore a different way is offered to solve this problem only using two thermo-infrared bands data based on three assumptions, constant emissivities in two measurements, and the same atmospheric parameters for neighbouring pixels and the difference of emissivity (Δε) of two channels can be known beforehand. Results of digital simulations show that it is possible to retrieve LST with its root mean square (RMS) of errors less than 1 K and RMS of relative error of ground radiance at 7% if the error of atmospheric temperature at ± 2°C and the relative error of atmospheric water vapor at ± 10% can be satisfied. Results have been confirmed by initial field test. Project supported by the National Natural Science Foundation of China (Grant No. 49471056) and China’s National Key Basic Research Plan.     

Correlation-based temperature and emissivity separation algorithm

Based on analyzing the relationship between the atmospheric downward radiance and surface emis- sivity, this paper proposes a correlation criterion to optimize surface temperature during the process of temperature and emissivity separation from thermal infrared hyperspectral data, and puts forward the correlation-based temperature and emissivity separation algorithm (CBTES). The algorithm uses the correlation between the atmospheric downward radiance and surface emissivity to optimize surface temperature, and obtains surface emissivity with this temperature. The accuracy of CBTES was evalu- ated by the simulated thermal infrared hyperspectral data. The simulated results show that the CBTES can achieve high accuracy of temperature and emissivity inversion. CBTES has been compared with the iterative spectrally smooth temperature/emissivity separation (ISSTES), and the comparison results show that they have relative accuracy. Besides, CBTES is insensitive to the instrumental random noise and the change of atmospheric downward radiance during the measurements. As regards the noniso- thermal pixel, its radiometric temperature changes slowly with the wavenumber when its emissivity is defined as r-emissivity. The CBTES can be used to derive the equivalent temperature of nonisothermal pixel in a narrow spectral region when we assumed that the radiometric temperature is invariable in the narrow spectral region. The derived equivalent temperatures in multi-spectral regions in 714―1250 cm?1 can characterize the change trend of nonisothermal pixel's radiometric temperature.     

月表有效太阳辐照度实时模型

月表太阳辐射是深入研究月表温度分布的关键问题之一. 本文根据月表有效太阳辐照度与太阳常数、太阳辐射入射角以及日月距离之间的关系,建立月表有效太阳辐照度的实时模型. 该模型在1950～2050年的100年内的理论误差百分比小于0.28％, 相对前人提出的模型在精度上有了很大程度的提高,能较为真实地反映月表有效太阳辐照度随时间的变化规律. 计算结果表明2007年月表太阳辐照度的年变化范围在1321.5～1416.6 W·m^-2之间,平均为1368.0 W·m^-2. 通过对月表太阳辐射入射角计算结果的分析,证实了月球的两极可能存在极昼极夜.     

Lunar Laser Ranging: Glorious Past And A Bright Future

Lunar Laser Ranging (LLR), a part of the NASA Apollo program, has beenon-going for more than 30 years. It provides the grist for a multi-disciplinarydata analysis mill. Results exist for solid Earth sciences, geodesy and geodynamics,solar system ephemerides, terrestrial and celestial reference frames, lunar physics,general relativity and gravitational theory. Combined with other data, it treatsprecession of the Earth's spin axis, lunar induced nutation, polar motion/Earthrotation, Earth orbit obliquity to the ecliptic, intersection of the celestial equatorwith the ecliptic, luni-solar solid body tides, lunar tidal deceleration, lunar physicaland free librations, structure of the moon and energy dissipation in the lunar interior.LLR provides input to lunar surface cartography and surveying, Earth station and lunar retroreflector location and motion, mass of the Earth-moon system, lunar and terrestrial gravity harmonics and Love numbers, relativistic geodesic precession, and the equivalence principle of general relativity. With the passive nature of the reflectors and steady improvement in observing equipment and data analysis, LLR continues to provide state-of-the-art results. Gains are steady as the data-base expands. After more than 30 years, LLR remains the only active Apollo experiment. It is important to recognize examples of efficient and cost effective progress of research. LLR is just such an example.     

On possible drivers of Sun-induced climate changes

We tested the validity of two current hypotheses on the dependence of climate change on solar activity. One of them states that variations in the tropospheric temperature are caused directly by changes of the solar radiance (total or spectral). The other suggests that cosmic ray (CR) fluctuations, caused by the solar/heliospheric modulation, affect the climate via cloud formation. Confronting these hypotheses with seven different sets of the global/hemispheric temperature reconstructions for the last 400 years, we found that the former mechanism is in general more prominent than the latter. Therefore, we can conclude that in so far as the Sun–climate connection is concerned tropospheric temperatures are more likely affected by variations in the UV radiation flux rather than by those in the CR flux.     

Remote measurement of river morphology via fusion of LiDAR topography and spectrally based bathymetry

This study developed and evaluated a hybrid approach to remote measurement of river morphology that combines LiDAR topography with spectrally based bathymetry. Comparison of filtered LiDAR point clouds with surveyed cross‐sections indicated that subtle features on low‐relief floodplains were accurately resolved by LiDAR but that submerged areas could not be detected due to strong absorption of near‐infrared laser pulses by water. The reduced number of returns made the active channel evident in a LiDAR point density map. A second dataset suggested that pulse intensity also could be used to discriminate land from water via a threshold‐based masking procedure. Fusion of LiDAR and optical data required accurate co‐registration of images to the LiDAR, and we developed an object‐oriented procedure for achieving this alignment. Information on flow depths was derived by correlating pixel values with field measurements of depth. Highly turbid conditions dictated a positive relation between green band radiance and flow depth and contributed to under‐prediction of pool depths. Water surface elevations extracted from the LiDAR along the water's edge were used to produce a continuous water surface that preserved along‐channel variations in slope. Subtracting local flow depths from this surface yielded estimates of the bed elevation that were then combined with LiDAR topography for exposed areas to create a composite representation of the riverine terrain. The accuracy of this terrain model was assessed via comparison with detailed field surveys. A map of elevation residuals showed that the greatest errors were associated with underestimation of pool depths and failure to capture cross‐stream differences in water surface elevation. Nevertheless, fusion of LiDAR and passive optical image data provided an efficient means of characterizing river morphology that would not have been possible if either dataset had been used in isolation. Copyright © 2011 John Wiley & Sons, Ltd.     

基于球冠域的月球重力地形校正方法

利用卫星重力测量资料研究地球、月球与火星等星体的内部构造时,需要进行重力地形校正,计算全球布格重力异常,而在球坐标中实现地形校正计算是一种有效的途径.本文提出球坐标系中的球冠域地形校正计算方法,给出了该方法涉及的球坐标系之间坐标转换方法和球冠域内地形模型重构方法,并进行理论验证．作者利用嫦娥一号激光测高数据对月球重力进行地形校正,获得了月球全球布格重力异常,并与球坐标系中Tesseroid 单元体地形校正方法对比,分析了球冠域地形校正方法的计算精度、空间分辨率及其优缺点．     

Thermal regime of a headwater stream within a clear‐cut,coastal British Columbia,Canada

This study examined the thermal regime of a headwater stream within a clear‐cut. The stream had a complex morphology dominated by step–pool features, many formed by sediment accumulation upstream of woody debris. Maximum daily temperatures increased up to 5 °C after logging, and were positively associated with maximum daily air temperature and negatively with discharge. Maximum daily temperatures generally increased with downstream distance through the cut block, but decreased with distance in two segments over distances of tens of metres, where the topography indicated relatively concentrated lateral inflow. Localized cool areas within a step–pool unit were associated with zones of concentrated upwelling. Bed temperatures tended to be higher and have greater ranges in areas of downwelling flow into the bed. Heat budget estimates were made using meteorological measurements over the water surface and a model of net radiation using canopy characteristics derived from fisheye photography. Heat exchange driven by hyporheic flow through the channel step was a cooling effect during daytime, with a magnitude up to approximately 25% that of net radiation during the period of maximum daytime warming. Heat budget calculations in these headwater streams are complicated by the heterogeneity of incident solar radiation and channel geometry, as well as uncertainty in estimating heat and water exchanges between the stream and the subsurface via hyporheic exchange and heat conduction. Copyright © 2005 John Wiley & Sons, Ltd.