Code and readme of a MATLAB-based graphical user interface program for computing functionals of the geopotential

We present a novel graphical user interface program GrafLab (GRAvity Field LABoratory) for spherical harmonic synthesis (SHS) created in MATLAB®. This program allows to comfortably compute 38 various functionals of the geopotential up to ultra-high degrees and orders of spherical harmonic expansion. For the most difficult part of the SHS, namely the evaluation of the fully normalized associated Legendre functions (fnALFs), we used three different approaches according to required maximum degree: (i) the standard forward column method (up to maximum degree 1800, in some cases up to degree 2190); (ii) the modified forward column method combined with Horner's scheme (up to maximum degree 2700); (iii) the extended-range arithmetic (up to an arbitrary maximum degree). For the maximum degree 2190, the SHS with fnALFs evaluated using the extended-range arithmetic approach takes only approximately 2-3 times longer than its standard arithmetic counterpart, i.e. the standard forward column method. In the GrafLab, the functionals of the geopotential can be evaluated on a regular grid or point-wise, while the input coordinates can either be read from a data file or entered manually. For the computation on a regular grid we decided to apply the lumped coefficients approach due to significant time-efficiency of this method. Furthermore, if a full variance-covariance matrix of spherical harmonic coefficients is available, it is possible to compute the commission errors of the functionals. When computing on a regular grid, the output functionals or their commission errors may be depicted on a map using automatically selected cartographic projection.

Modelling calving front dynamics using a level-set method: application to Jakobshavn Isbræ, West Greenland, links model results in MATLAB format

Calving is a major mechanism of ice discharge of the Antarctic and Greenland ice sheets, and a change in calving front position affects the entire stress regime of marine terminating glaciers. The representation of calving front dynamics in a 2-D or 3-D ice sheet model remains non-trivial. Here, we present the theoretical and technical framework for a level-set method, an implicit boundary tracking scheme, which we implement into the Ice Sheet System Model (ISSM). This scheme allows us to study the dynamic response of a drainage basin to user-defined calving rates. We apply the method to Jakobshavn Isbræ, a major marine terminating outlet glacier of the West Greenland Ice Sheet. The model robustly reproduces the high sensitivity of the glacier to calving, and we find that enhanced calving triggers significant acceleration of the ice stream. Upstream acceleration is sustained through a combination of mechanisms. However, both lateral stress and ice influx stabilize the ice stream. This study provides new insights into the ongoing changes occurring at Jakobshavn Isbræ and emphasizes that the incorporation of moving boundaries and dynamic lateral effects, not captured in flow-line models, is key for realistic model projections of sea level rise on centennial timescales.

Set of Matlab function to study geometrical parameters or channels

This set of functions allows one to compute the radius of curvature of a river in planform for the purpose of making correlations with other geometric parameters of a channel. The code may also be used to compute the width of a channel.

Matlab scripts for modeled obliquity amplitude modulations

The start of the Mesozoic Era is marked by roughly five million years (myr) of Earth system upheavals, including unstable biotic recovery, repeated global warming, ocean anoxia, and perturbations in the global carbon cycle. Intervals between crises were comparably hospitable to life. The causes of these upheavals are unknown, but are thought to be linked to recurrent Siberian volcanism. Here, two marine sedimentary successions at Chaohu and Daxiakou, South China are evaluated for paleoclimate change from astronomical forcing. In these sections, gamma-ray variations indicative of terrestrial weathering reveal enhanced obliquity cycling over prolonged intervals, characterized by a periodicity of 32.8 kiloyear and strong 1.2 myr modulations. This suggests a 22-hour length-of-day and 1.2 myr interaction between the orbital inclinations of Earth and Mars. The 1.2 myr obliquity modulation cycles in these sections are compared with Early Triassic records of global sea-level, temperature, redox and biotic evolution. The evidence collectively suggests that long-term astronomical forcing was involved in the repeated climatic and biotic upheavals that took place throughout the Early Triassic.

Directional wave measurements using the Offshore Sensing Sailbuoy, link to MATLAB file

An autonomous vessel, the Offshore Sensing Sailbuoy, was used for wave measurements near the Ekofisk oil platform complex in the North Sea (56.5 N, 3.2 E, operated by ConocoPhilllips) from 6 to 20 November 2015. Being 100% wind propelled, the Sailbuoy has two-way communication via the Iridium network and has the capability for missions of six months or more. It has previously been deployed in the Arctic, Norwegian Sea and the Gulf of Mexico, but this was the first real test for wave measurements. During the campaign it held position about 20km northeast of Ekofisk (on the lee side) during rough conditions. Mean wind speed measured at Ekofisk during the campaign was near 9.8m/s, with a maximum of 20.4m/s, with wind mostly from south and south west. A Datawell MOSE G1000 GPS based 2Hz wave sensor was mounted on the Sailbuoy. Mean significant wave height (Hs 1hr) measured was 3m, whereas maximum Hs was 6m. Mean wave period was 7.7s, while maximum wave height, Hmax, was 12.6m. These measurements have been compared with non-directional Waverider observations at the Ekofisk complex. Mean Hs at Ekofisk was 3.1m, while maximum Hs was 6.5m. Nevertheless, the correlation between the two measurements was high (97%). Spectra comparison was also good, except for low Hs (~1m), where the motion of the vessel seemed to influence the measurements. Nevertheless, the Sailbuoy performed well during this campaign, and results suggests that it is a suitable platform for wave measurements in rather rough sea conditions.

Potential temperature in Espartel Sill (Strait of Gibraltar), from 2004-09-30 to 2016-03-02 at 344 meters depth, link to raw data in MATLAB format

Potential temperature measured with a SBE37 at 35.862ºN, 5.97ºW at 344 meters Depth. Data expand from September the 30th, 2004 to March the 2nd, 2016. Original measurement frequency was 30 minutes, the data presented here is a subsampling that extract the coldest temperature found each 12 hours. The time vector corresponds with the moment in which this minimun temperature is observed.