Analemma ও Solar Declination: আকাশে সূর্যের বার্ষিক ‘৮’ আকৃতি কেন?

Analemma একটি সুন্দর photograph মাত্র নয়; এটি একই diagram-এ Earth-এর axial tilt, orbital eccentricity, solar declination এবং apparent solar time-এর অনিয়মকে দৃশ্যমান করে। Panchanga software-এর জন্য এটি একটি শক্তিশালী diagnostic—কারণ solar position-এর north–south এবং east–west behavior একই সঙ্গে পরীক্ষা করা যায়।

Analemma কী?

ধরা যাক, একটি স্থির camera দিয়ে প্রতি সপ্তাহে দুপুর 12:00-এ সূর্যের ছবি তোলা হলো। Camera, lens, frame, direction এবং clock convention অপরিবর্তিত রাখা হলো। এক বছর পরে ছবিগুলো একই frame-এ composite করলে Sun-এর বিন্দুগুলো একটি ‘৮’-এর মতো curve তৈরি করবে।

NASA-র Astronomy Picture of the Day-এর বহু composite-এ ঠিক এই পদ্ধতি দেখা যায়—একই স্থান থেকে একই clock time-এ নেওয়া বহু exposure একত্র করা হয়। Solstice-এর Sun curve-এর উত্তরতম ও দক্ষিণতম অংশে থাকে; equinox-এর Sun মাঝামাঝি আসে, যদিও সাধারণত figure-eight-এর crossing point-এ নয়।

এটি Sun-এর orbit নয়

Analemma-কে আকাশে Sun-এর বাস্তব orbital path ভাবা ভুল। Earth থেকে Sun-এর বার্ষিক background-star path হলো ecliptic। Analemma হলো একটি বিশেষ sampling pattern: প্রতি solar day-এর একই clock time-এ Sun-এর position নেওয়া। Sampling rule বদলালে curve-ও বদলে যায়।

Solar Declination কী?

Declination হলো celestial equator-এর উত্তর বা দক্ষিণে কোনো celestial body-এর angular distance—পৃথিবীর latitude-এর আকাশীয় সমতুল্য। Sun-এর geocentric declination সাধারণত δ দ্বারা লেখা হয়:

• δ ধনাত্মক হলে Sun celestial equator-এর উত্তরে;
• δ ঋণাত্মক হলে Sun celestial equator-এর দক্ষিণে;
• δ = 0°-এর কাছে Sun celestial equator অতিক্রম করে;
• বর্তমান যুগে extrema প্রায় +23.4° ও −23.4°।

এই সীমা Earth-এর obliquity বা axial tilt-এর কারণে। June solstice-এর কাছে declination উত্তরতম, December solstice-এর কাছে দক্ষিণতম এবং March ও September equinox-এর কাছে শূন্য হয়।

‘৮’ আকৃতির দুই অক্ষ

Analemma-র shape বোঝার সবচেয়ে সহজ উপায় হলো curve-টিকে দুই independent-looking component-এ ভাগ করা:

উল্লম্ব বা উত্তর–দক্ষিণ component

এটি মূলত solar declination-এর seasonal variation। Northern summer-এ Sun আকাশে উত্তরে ও উঁচুতে, winter-এ দক্ষিণে ও নিচুতে থাকে। Camera frame-এ vertical direction exact celestial north–south না-ও হতে পারে, কিন্তু seasonal height shift-এর প্রধান কারণ declination।

অনুভূমিক বা পূর্ব–পশ্চিম component

একই clock time-এ true Sun কখনও mean Sun-এর আগে, কখনও পরে থাকে। এই apparent-minus-mean time difference হলো Equation of Time। আকাশে সময়ের ১ মিনিট প্রায় 0.25° hour angle-এর সমতুল্য; তাই ±16 মিনিটের variation প্রায় ±4° east–west displacement দিতে পারে।

বাস্তব photograph-এ horizon coordinate projection, camera direction, latitude এবং chosen clock time এই দুই component-কে rotate, shear বা distort করতে পারে। তাই সব analemma vertical ‘৮’ হিসেবে দেখা যাবে না।

দুটি loop সমান নয় কেন?

Earth-এর orbit যদি perfectly circular হতো, orbital angular speed আরও uniform হতো। আবার axis যদি tilted না হতো, declination variation থাকত না। বাস্তবে axial tilt ও orbital eccentricity দুটিই আছে এবং তাদের phase একই নয়।

1. Obliquity effect: ecliptic motion-এর projection celestial equator-এ uniform নয়;
2. Eccentricity effect: perihelion-এর কাছে Earth দ্রুত, aphelion-এর কাছে ধীরে চলে;
3. Phase relation: দুই effect বছরের আলাদা অংশে maxima/minima পায়;
4. Combined result: east–west displacement এবং north–south displacement symmetric হয় না।

ফলে analemma-র একটি loop বড় এবং অন্যটি ছোট হয়। Shape দীর্ঘ সময়ে অপরিবর্তনীয়ও নয়—obliquity, eccentricity এবং perihelion orientation ধীরে বদলায়। Historical millennia-র জন্য আধুনিক photograph-এর shape হুবহু ধরে নেওয়া যায় না।

Solstice, Equinox ও Crossing Point

Solstice-এ declination extrema হওয়ায় analemma-র উত্তরতম ও দক্ষিণতম point পাওয়া যায়। কিন্তু equinox-এ declination zero হলেও analemma-র crossing point হওয়া আবশ্যক নয়।

Figure-eight crossing-এর জন্য দুই ভিন্ন date-এ simultaneously একই declination এবং একই Equation of Time দরকার। Equinox শুধু declination zero হওয়ার কাছাকাছি seasonal condition জানায়; Equation of Time তখন zero হবে—এমন নিয়ম নেই। তাই September equinox-এর Sun analemma-র মাঝামাঝি থাকলেও crossing point থেকে সরে থাকতে পারে।

Analemma Photograph কীভাবে তৈরি হয়?

একটি নির্ভরযোগ্য photographic analemma তৈরিতে astronomy-এর পাশাপাশি কঠোর observational consistency দরকার:

• একই camera position ও fixed tripod mount;
• একই focal length, framing ও orientation;
• একই clock time এবং ঘোষিত DST policy;
• নিরাপদ certified solar filter;
• নিয়মিত interval—যেমন প্রতি ৭ বা ১০ দিনে;
• এক বছরের date/time log;
• একটি পৃথক foreground exposure;
• cloudy day ও missing exposure-এর স্বচ্ছ record।

Fixed clock time বনাম fixed apparent solar time

DST চালু হওয়ার দিনে clock time এক ঘণ্টা বদলে গেলে curve-এ discontinuity তৈরি হতে পারে। Scientific analemma-য় সাধারণত year-round standard time, fixed UTC অথবা explicitly corrected local mean time ব্যবহার করা ভালো।

Solar Noon-এ Sun কত উঁচুতে?

Geometric meridian altitude-এর একটি সহজ relation:

এই সূত্র upper culmination-এর geometric altitude দেয় এবং refraction, observer elevation ও সূর্যের disk size উপেক্ষা করে। δ = φ হলে Sun zenith দিয়ে transit করে। Tropical latitude-এর মধ্যে বছরে সাধারণত দুবার এমন হতে পারে; tropic-এর ঠিক বাইরে কখনও exact zenith transit হয় না।

Sun কি সব সময় দক্ষিণে transit করে?

না। Northern tropics-এর কোনো observer-এর latitude-এর চেয়ে solar declination বেশি উত্তরে হলে noon Sun zenith-এর উত্তর দিক দিয়ে transit করে। তাই “দুপুরে Sun সব সময় দক্ষিণে”—এটি mid-latitude intuition, universal rule নয়।

Declination ও দিনের দৈর্ঘ্য

Ideal point-Sun ও zero-altitude horizon model-এ daylight duration-এর approximation:

এখানে inverse cosine degree-এ নিতে হবে। বাস্তব sunrise/sunset-এর জন্য solar upper limb, atmospheric refraction, elevation ও terrain horizon যোগ করতে হয়। Polar circle-এর বাইরে arccos argument range-এর বাইরে গেলে ordinary rise/set নেই—সেটিকে error নয়, physical status হিসেবে handle করতে হবে।

Declination একই সঙ্গে sunrise azimuth, sunset azimuth, day length, noon altitude এবং twilight duration বদলায়। Equation of Time day length তৈরি করে না; এটি sunrise, transit ও sunset-এর clock label-কে একসঙ্গে আগে বা পরে সরায়।

কলকাতা, ঢাকা ও নিউইয়র্কে seasonal altitude

নিচের values geometric এবং rounded; refraction ধরা হয়নি।

কলকাতা Tropic of Cancer-এর দক্ষিণে, তাই Sun বছরে দুইবার zenith-এর খুব কাছ দিয়ে যেতে পারে। ঢাকা বর্তমান tropic-এর সামান্য উত্তরে; তাই June solstice-এ Sun অত্যন্ত উঁচু হলেও exact zenith transit সাধারণত হয় না। নিউইয়র্কে seasonal altitude range অনেক বেশি দৃশ্যমান এবং winter noon Sun তুলনামূলক নিচু।

Panchanga Calculation-এ এর ব্যবহার

Solar declination নিজে বাংলা সৌর মাসের boundary নয়, কিন্তু পঞ্জিকার বহু স্থানীয় ঘটনার geometry-তে অপরিহার্য:

• sunrise ও sunset hour angle;
• দিনমান ও রাত্রিমান;
• civil, nautical ও astronomical twilight;
• local solar-noon altitude ও shadow direction;
• Arunodaya ও daytime segment-এর clock time;
• high-latitude no-rise/no-set detection;
• observational validation ও annual solar charts।

Declination আর sidereal longitude এক নয়

Sun-এর declination zero দেখে ১ বৈশাখ বা Mesha Saṅkrānti ঘোষণা করা যাবে না। বাংলা সৌর মাসের জন্য project-এর ঘোষিত nirayana/Surya Siddhanta longitude ও month-start rule ব্যবহার করতে হবে।

C#-এ Approximate Analemma Data

NOAA-এর compact formula দিয়ে educational chart বা sanity check তৈরি করা যায়। High-precision production result-এর জন্য ephemeris থেকে apparent right ascension ও declination নিন।

static double FractionalYear(DateTime date, double localHour)
{
    int days = DateTime.IsLeapYear(date.Year) ? 366 : 365;
    return 2.0 * Math.PI / days *
        (date.DayOfYear - 1 + (localHour - 12.0) / 24.0);
}

static double SolarDeclinationRadians(double gamma)
{
    return 0.006918
         - 0.399912 * Math.Cos(gamma)
         + 0.070257 * Math.Sin(gamma)
         - 0.006758 * Math.Cos(2 * gamma)
         + 0.000907 * Math.Sin(2 * gamma)
         - 0.002697 * Math.Cos(3 * gamma)
         + 0.001480 * Math.Sin(3 * gamma);
}

static double EquationOfTimeMinutes(double gamma)
{
    return 229.18 *
        (0.000075
        + 0.001868 * Math.Cos(gamma)
        - 0.032077 * Math.Sin(gamma)
        - 0.014615 * Math.Cos(2 * gamma)
        - 0.040849 * Math.Sin(2 * gamma));
}

Analemma point model

public sealed class AnalemmaPoint
{
    public DateTime Date { get; set; }
    public double DeclinationDegrees { get; set; }
    public double EquationOfTimeMinutes { get; set; }
}

static AnalemmaPoint BuildPoint(DateTime date, double clockHour)
{
    double gamma = FractionalYear(date, clockHour);
    return new AnalemmaPoint
    {
        Date = date.Date,
        DeclinationDegrees =
            SolarDeclinationRadians(gamma) * 180.0 / Math.PI,
        EquationOfTimeMinutes = EquationOfTimeMinutes(gamma)
    };
}

Simple scientific plot-এ x-axis হিসেবে Equation of Time এবং y-axis হিসেবে declination ব্যবহার করলে location-independent diagnostic analemma পাওয়া যায়। Photographic horizon plot তৈরি করতে observer latitude, longitude, time zone, hour angle, refraction policy এবং azimuth/altitude projection দরকার।

Noon altitude function

static double GeometricNoonAltitude(
    double latitudeDegrees,
    double declinationDegrees)
{
    return 90.0 - Math.Abs(
        latitudeDegrees - declinationDegrees);
}

Function-টি geometric upper culmination-এর জন্য। Apparent altitude দেখাতে refraction যোগ করতে হবে; lower culmination বা polar cases আলাদা status প্রয়োজন।

Validation Checklist

• একই date/time scale-এ declination ও E গণনা করা হয়েছে;
• Degrees ও radians মিশে যায়নি;
• Leap year denominator 366 ব্যবহার করা হয়েছে;
• Equation of Time sign প্রকাশ্যে documented;
• June solstice-এ declination প্রায় +23.4°;
• December solstice-এ declination প্রায় −23.4°;
• Equinox-এর কাছে declination zero crossing test করা হয়েছে;
• Analemma crossing-কে equinox ধরে নেওয়া হয়নি;
• Fixed clock time এবং fixed solar-noon sampling আলাদা;
• DST discontinuity chart metadata-তে উল্লেখ করা হয়েছে;
• Solar declination ও sidereal longitude পৃথক field;
• Approximate formula ephemeris values-এর সঙ্গে তুলনা করা হয়েছে;
• Polar invalid arccos typed status হিসেবে ফেরে;
• Historical year-এ obliquity/date model উল্লেখ করা হয়েছে।

সাধারণ ভুল

1. Analemma-কে Sun-এর orbit বলা: এটি fixed-time sampling pattern।
2. Equinox-কে crossing point বলা: E zero না হলে crossing হয় না।
3. Clock noon ও solar noon মেশানো: Equation of Time component হারিয়ে যায়।
4. Declination দিয়ে Saṅkrānti নির্ণয়: ecliptic longitude-এর পরিবর্তে equatorial coordinate ব্যবহৃত হয়।
5. DST-সহ raw local clock ব্যবহার: curve-এ artificial one-hour jump আসে।
6. Altitude-কে declination বলা: altitude observer latitude ও hour angle-এর উপরও নির্ভরশীল।
7. Approximation-কে ephemeris precision বলা: educational formula-এর সীমা গোপন থাকে।

উপসংহার

Analemma হলো Earth–Sun geometry-এর একটি বার্ষিক স্বাক্ষর। Solar declination সূর্যের north–south seasonal migration তৈরি করে; Equation of Time fixed clock-এর তুলনায় east–west displacement তৈরি করে। Axial tilt ও elliptical orbit-এর যৌথ প্রভাবে দুটি component মিলিয়ে অসম figure-eight গঠিত হয়।

Panchanga-তে analemma সরাসরি কোনো তিথি বা মাস নির্ধারণ করে না, কিন্তু solar model যাচাইয়ের অসাধারণ উপায়। একই engine-এর declination, Equation of Time, solar noon, sunrise–sunset ও daylight duration পরস্পরের সঙ্গে সঙ্গতিপূর্ণ কি না—একটি annual analemma report তা দ্রুত প্রকাশ করতে পারে।

তথ্যসূত্র ও আরও পাঠ

1. U.S. Naval Observatory—Sunrise and Sunset Times Near the Solstices; declination range, daylight formula, Equation of Time এবং analemma।
2. U.S. Naval Observatory—The Equation of Time; apparent/mean solar time, obliquity ও eccentricity effect।
3. NOAA—General Solar Position Calculations; approximate declination ও Equation of Time equations।
4. NASA Science—An Analemma of the Sun; same-time composite ও solstice positions-এর photographic example।
5. U.S. Naval Observatory—Earth’s Seasons; equinox, solstice, perihelion ও aphelion data service।