\documentclass[11pt,twoside]{article}\makeatletter

\IfFileExists{xcolor.sty}%
  {\RequirePackage{xcolor}}%
  {\RequirePackage{color}}
\usepackage{colortbl}
\usepackage{wrapfig}
\usepackage{ifxetex}
\ifxetex
  \usepackage{fontspec}
  \usepackage{xunicode}
  \catcode`⃥=\active \def⃥{\textbackslash}
  \catcode`❴=\active \def❴{\{}
  \catcode`❵=\active \def❵{\}}
  \def\textJapanese{\fontspec{Noto Sans CJK JP}}
  \def\textChinese{\fontspec{Noto Sans CJK SC}}
  \def\textKorean{\fontspec{Noto Sans CJK KR}}
  \setmonofont{DejaVu Sans Mono}
  
\else
  \IfFileExists{utf8x.def}%
   {\usepackage[utf8x]{inputenc}
      \PrerenderUnicode{–}
    }%
   {\usepackage[utf8]{inputenc}}
  \usepackage[english]{babel}
  \usepackage[T1]{fontenc}
  \usepackage{float}
  \usepackage[]{ucs}
  \uc@dclc{8421}{default}{\textbackslash }
  \uc@dclc{10100}{default}{\{}
  \uc@dclc{10101}{default}{\}}
  \uc@dclc{8491}{default}{\AA{}}
  \uc@dclc{8239}{default}{\,}
  \uc@dclc{20154}{default}{ }
  \uc@dclc{10148}{default}{>}
  \def\textschwa{\rotatebox{-90}{e}}
  \def\textJapanese{}
  \def\textChinese{}
  \IfFileExists{tipa.sty}{\usepackage{tipa}}{}
\fi
\def\exampleFont{\ttfamily\small}
\DeclareTextSymbol{\textpi}{OML}{25}
\usepackage{relsize}
\RequirePackage{array}
\def\@testpach{\@chclass
 \ifnum \@lastchclass=6 \@ne \@chnum \@ne \else
  \ifnum \@lastchclass=7 5 \else
   \ifnum \@lastchclass=8 \tw@ \else
    \ifnum \@lastchclass=9 \thr@@
   \else \z@
   \ifnum \@lastchclass = 10 \else
   \edef\@nextchar{\expandafter\string\@nextchar}%
   \@chnum
   \if \@nextchar c\z@ \else
    \if \@nextchar l\@ne \else
     \if \@nextchar r\tw@ \else
   \z@ \@chclass
   \if\@nextchar |\@ne \else
    \if \@nextchar !6 \else
     \if \@nextchar @7 \else
      \if \@nextchar (8 \else
       \if \@nextchar )9 \else
  10
  \@chnum
  \if \@nextchar m\thr@@\else
   \if \@nextchar p4 \else
    \if \@nextchar b5 \else
   \z@ \@chclass \z@ \@preamerr \z@ \fi \fi \fi \fi
   \fi \fi  \fi  \fi  \fi  \fi  \fi \fi \fi \fi \fi \fi}
\gdef\arraybackslash{\let\\=\@arraycr}
\def\@textsubscript#1{{\m@th\ensuremath{_{\mbox{\fontsize\sf@size\z@#1}}}}}
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\def\abbr{}
\def\corr{}
\def\expan{}
\def\gap{}
\def\orig{}
\def\reg{}
\def\ref{}
\def\sic{}
\def\persName{}\def\name{}
\def\placeName{}
\def\orgName{}
\def\textcal#1{{\fontspec{Lucida Calligraphy}#1}}
\def\textgothic#1{{\fontspec{Lucida Blackletter}#1}}
\def\textlarge#1{{\large #1}}
\def\textoverbar#1{\ensuremath{\overline{#1}}}
\def\textquoted#1{‘#1’}
\def\textsmall#1{{\small #1}}
\def\textsubscript#1{\@textsubscript{\selectfont#1}}
\def\textxi{\ensuremath{\xi}}
\def\titlem{\itshape}
\newenvironment{biblfree}{}{\ifvmode\par\fi }
\newenvironment{bibl}{}{}
\newenvironment{byline}{\vskip6pt\itshape\fontsize{16pt}{18pt}\selectfont}{\par }
\newenvironment{citbibl}{}{\ifvmode\par\fi }
\newenvironment{docAuthor}{\ifvmode\vskip4pt\fontsize{16pt}{18pt}\selectfont\fi\itshape}{\ifvmode\par\fi }
\newenvironment{docDate}{}{\ifvmode\par\fi }
\newenvironment{docImprint}{\vskip 6pt}{\ifvmode\par\fi }
\newenvironment{docTitle}{\vskip6pt\bfseries\fontsize{22pt}{25pt}\selectfont}{\par }
\newenvironment{msHead}{\vskip 6pt}{\par}
\newenvironment{msItem}{\vskip 6pt}{\par}
\newenvironment{rubric}{}{}
\newenvironment{titlePart}{}{\par }

\newcolumntype{L}[1]{){\raggedright\arraybackslash}p{#1}}
\newcolumntype{C}[1]{){\centering\arraybackslash}p{#1}}
\newcolumntype{R}[1]{){\raggedleft\arraybackslash}p{#1}}
\newcolumntype{P}[1]{){\arraybackslash}p{#1}}
\newcolumntype{B}[1]{){\arraybackslash}b{#1}}
\newcolumntype{M}[1]{){\arraybackslash}m{#1}}
\definecolor{label}{gray}{0.75}
\def\unusedattribute#1{\sout{\textcolor{label}{#1}}}
\DeclareRobustCommand*{\xref}{\hyper@normalise\xref@}
\def\xref@#1#2{\hyper@linkurl{#2}{#1}}
\begingroup
\catcode`\_=\active
\gdef_#1{\ensuremath{\sb{\mathrm{#1}}}}
\endgroup
\mathcode`\_=\string"8000
\catcode`\_=12\relax

\usepackage[a4paper,twoside,lmargin=1in,rmargin=1in,tmargin=1in,bmargin=1in,marginparwidth=0.75in]{geometry}
\usepackage{framed}

\definecolor{shadecolor}{gray}{0.95}
\usepackage{longtable}
\usepackage[normalem]{ulem}
\usepackage{fancyvrb}
\usepackage{fancyhdr}
\usepackage{graphicx}
\usepackage{marginnote}

\renewcommand{\@cite}[1]{#1}


\renewcommand*{\marginfont}{\itshape\footnotesize}

\def\Gin@extensions{.pdf,.png,.jpg,.mps,.tif}

  \pagestyle{fancy}

\usepackage[pdftitle={Scan to know paper details and author's profile 465U},
 pdfauthor={}]{hyperref}
\hyperbaseurl{}

	 \paperwidth210mm
	 \paperheight297mm
              
\def\@pnumwidth{1.55em}
\def\@tocrmarg {2.55em}
\def\@dotsep{4.5}
\setcounter{tocdepth}{3}
\clubpenalty=8000
\emergencystretch 3em
\hbadness=4000
\hyphenpenalty=400
\pretolerance=750
\tolerance=2000
\vbadness=4000
\widowpenalty=10000

\renewcommand\section{\@startsection {section}{1}{\z@}%
     {-1.75ex \@plus -0.5ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large\bfseries}}
\renewcommand\subsection{\@startsection{subsection}{2}{\z@}%
     {-1.75ex\@plus -0.5ex \@minus- .2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large}}
\renewcommand\subsubsection{\@startsection{subsubsection}{3}{\z@}%
     {-1.5ex\@plus -0.35ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\large}}
\renewcommand\paragraph{\@startsection{paragraph}{4}{\z@}%
     {-1ex \@plus-0.35ex \@minus -0.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\normalsize}}
\renewcommand\subparagraph{\@startsection{subparagraph}{5}{\parindent}%
     {1.5ex \@plus1ex \@minus .2ex}%
     {-1em}%
     {\reset@font\normalsize\bfseries}}


\def\l@section#1#2{\addpenalty{\@secpenalty} \addvspace{1.0em plus 1pt}
 \@tempdima 1.5em \begingroup
 \parindent \z@ \rightskip \@pnumwidth 
 \parfillskip -\@pnumwidth 
 \bfseries \leavevmode #1\hfil \hbox to\@pnumwidth{\hss #2}\par
 \endgroup}
\def\l@subsection{\@dottedtocline{2}{1.5em}{2.3em}}
\def\l@subsubsection{\@dottedtocline{3}{3.8em}{3.2em}}
\def\l@paragraph{\@dottedtocline{4}{7.0em}{4.1em}}
\def\l@subparagraph{\@dottedtocline{5}{10em}{5em}}
\@ifundefined{c@section}{\newcounter{section}}{}
\@ifundefined{c@chapter}{\newcounter{chapter}}{}
\newif\if@mainmatter 
\@mainmattertrue
\def\chaptername{Chapter}
\def\frontmatter{%
  \pagenumbering{roman}
  \def\thechapter{\@roman\c@chapter}
  \def\theHchapter{\roman{chapter}}
  \def\thesection{\@roman\c@section}
  \def\theHsection{\roman{section}}
  \def\@chapapp{}%
}
\def\mainmatter{%
  \cleardoublepage
  \def\thechapter{\@arabic\c@chapter}
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \pagenumbering{arabic}
  \setcounter{secnumdepth}{6}
  \def\@chapapp{\chaptername}%
  \def\theHchapter{\arabic{chapter}}
  \def\thesection{\@arabic\c@section}
  \def\theHsection{\arabic{section}}
}
\def\backmatter{%
  \cleardoublepage
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \setcounter{secnumdepth}{2}
  \def\@chapapp{\appendixname}%
  \def\thechapter{\@Alph\c@chapter}
  \def\theHchapter{\Alph{chapter}}
  \appendix
}
\newenvironment{bibitemlist}[1]{%
   \list{\@biblabel{\@arabic\c@enumiv}}%
       {\settowidth\labelwidth{\@biblabel{#1}}%
        \leftmargin\labelwidth
        \advance\leftmargin\labelsep
        \@openbib@code
        \usecounter{enumiv}%
        \let\p@enumiv\@empty
        \renewcommand\theenumiv{\@arabic\c@enumiv}%
	}%
  \sloppy
  \clubpenalty4000
  \@clubpenalty \clubpenalty
  \widowpenalty4000%
  \sfcode`\.\@m}%
  {\def\@noitemerr
    {\@latex@warning{Empty `bibitemlist' environment}}%
    \endlist}

\def\tableofcontents{\section*{\contentsname}\@starttoc{toc}}
\parskip0pt
\parindent1em
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\newenvironment{reflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemsep}{0pt}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\itshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{sansreflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\itemsep}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\upshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{specHead}[2]%
 {\vspace{20pt}\hrule\vspace{10pt}%
  \phantomsection\label{#1}\markright{#2}%

  \pdfbookmark[2]{#2}{#1}%
  \hspace{-0.75in}{\bfseries\fontsize{16pt}{18pt}\selectfont#2}%
  }{}
      \def\TheFullDate{1970-01-01 (revised: 01 January 1970)}
\def\TheID{\makeatother }
\def\TheDate{1970-01-01}
\title{Scan to know paper details and author's profile 465U}
\author{}\makeatletter 
\makeatletter
\newcommand*{\cleartoleftpage}{%
  \clearpage
    \if@twoside
    \ifodd\c@page
      \hbox{}\newpage
      \if@twocolumn
        \hbox{}\newpage
      \fi
    \fi
  \fi
}
\makeatother
\makeatletter
\thispagestyle{empty}
\markright{\@title}\markboth{\@title}{\@author}
\renewcommand\small{\@setfontsize\small{9pt}{11pt}\abovedisplayskip 8.5\p@ plus3\p@ minus4\p@
\belowdisplayskip \abovedisplayskip
\abovedisplayshortskip \z@ plus2\p@
\belowdisplayshortskip 4\p@ plus2\p@ minus2\p@
\def\@listi{\leftmargin\leftmargini
               \topsep 2\p@ plus1\p@ minus1\p@
               \parsep 2\p@ plus\p@ minus\p@
               \itemsep 1pt}
}
\makeatother
\fvset{frame=single,numberblanklines=false,xleftmargin=5mm,xrightmargin=5mm}
\fancyhf{} 
\setlength{\headheight}{14pt}
\fancyhead[LE]{\bfseries\leftmark} 
\fancyhead[RO]{\bfseries\rightmark} 
\fancyfoot[RO]{}
\fancyfoot[CO]{\thepage}
\fancyfoot[LO]{\TheID}
\fancyfoot[LE]{}
\fancyfoot[CE]{\thepage}
\fancyfoot[RE]{\TheID}
\hypersetup{citebordercolor=0.75 0.75 0.75,linkbordercolor=0.75 0.75 0.75,urlbordercolor=0.75 0.75 0.75,bookmarksnumbered=true}
\fancypagestyle{plain}{\fancyhead{}\renewcommand{\headrulewidth}{0pt}}

\date{}
\usepackage{authblk}

\providecommand{\keywords}[1]
{
\footnotesize
  \textbf{\textit{Index terms---}} #1
}

\usepackage{graphicx,xcolor}
\definecolor{GJBlue}{HTML}{273B81}
\definecolor{GJLightBlue}{HTML}{0A9DD9}
\definecolor{GJMediumGrey}{HTML}{6D6E70}
\definecolor{GJLightGrey}{HTML}{929497} 

\renewenvironment{abstract}{%
   \setlength{\parindent}{0pt}\raggedright
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
   \textcolor{GJBlue}{\large\bfseries\abstractname\space}
}{%   
   \vskip8pt
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
}

\usepackage[absolute,overlay]{textpos}

\makeatother 
      \usepackage{lineno}
      \linenumbers
      
\begin{document}

             \affil[1]{  }

\renewcommand\Authands{ and }

\date{\small \em Received: 1 January 1970 Accepted: 1 January 1970 Published: 1 January 1970}

\maketitle


\begin{abstract}
        




\end{abstract}


\keywords{}

\begin{textblock*}{18cm}(1cm,1cm) % {block width} (coords) 
\textcolor{GJBlue}{\LARGE Global Journals \LaTeX\ JournalKaleidoscope\texttrademark}
\end{textblock*}

\begin{textblock*}{18cm}(1.4cm,1.5cm) % {block width} (coords) 
\textcolor{GJBlue}{\footnotesize \\ Artificial Intelligence formulated this projection for compatibility purposes from the original article published at Global Journals. However, this technology is currently in beta. \emph{Therefore, kindly ignore odd layouts, missed formulae, text, tables, or figures.}}
\end{textblock*}


\begin{textblock*}{10cm}(1.05cm,3cm)
{{\textit{CrossRef DOI of original article:}} \underline{}}
\end{textblock*}\let\tabcellsep& 	 	 		 
\section[{I. INTRODUCCIÃ?"N}]{I. INTRODUCCIÃ?"N}\par
La hidrología y los SIG (Sistemas de información geográfica) se han constituidos como dos áreas que comparten una gran variedad de información \hyperref[b13]{(Morad, 2001)}. Por tal motivo, ha sido de gran interés el estudio de las relaciones y vínculos que existen entre los SIG y el modelamiento hidrológico \hyperref[b2]{( Castro López, 2011)}. El método del número de curva (CN) es uno de los vínculos que presentan estos dos campos de investigación. Esta metodología fue desarrollada por SCS (Soil Conservation Service) de los Estados Unidos con la finalidad de determinar la precipitación efectiva de una cuenca a partir del tipo de cobertura y uso de suelo, la lluvia acumulada y las condiciones de humedad de la zona \hyperref[b22]{(Ven Te, 1984)}; es decir, a partir de estas características físicas del territorio, este parámetro hidrológico permite establecer el potencial de escorrentía de la cuenca \hyperref[b7]{( Díaz Carvajal , 2017)}.\par
De acuerdo con \hyperref[b17]{(Romano Castillo, 2007)} el CN es el método más adecuado para transformar la precipitación total en precipitación efectiva debido a que tiene en cuenta una de las principales etapas del ciclo del agua, como lo es la infiltración, la cual influye notablemente en la dinámica hidrológica. Este proceso se presenta cuando el agua precipitada se mueve de manera vertical a través de la superficie del suelo generando el incremento en la humedad en esta matriz \hyperref[b1]{(Calvache, 2010)}.\par
A pesar de que la metodología del SCS se ha convertido en una de las principales herramientas para la modelación hidrológica y ambiental (Sánchez \hyperref[b19]{Palomino, 2015)}, el número de curva ha generado gran interés en las diferentes investigaciones debido a que permite establecer y describir las condiciones geográficas generales de una zona de estudio (Sedano Cruz, 2013) al tiempo que indica el potencial de escorrentía de una superficie determinada.\par
Para la estimación del número de curva en una cuenca hidrográfica, se han aplicado los SIG y la teledetección como herramientas de gran utilidad para determinar las condiciones generales de la cobertura vegetal y del uso del suelo(Lao \hyperref[b9]{Ramos, 2018)}.  {\ref (Deshmukh et.al., 2013)} emplearon y procesaron imágenes satelitales por medio de un software SIG con el fin de identificar la transformación del CN que se ha presentado en tres cuencas hidrográficas ubicadas en India. Por medio de los resultados obtenidos se pudo concluir que este parámetro hidrológico se ha incrementado en la zona a raíz del crecimiento desmesurado de la frontera agrícola.\par
Por su parte,  {\ref (Rodríguez et.al., 2012)} establecieron mediante la zonificación del CN, el impacto ambiental que se presenta en un escenario determinado, en cuanto a la ocupación del territorio, por medio de uso de las herramientas SIG. En Colombia, se llevó a cabo un estudio en la subcuenca de Betancí, ubicada en el departamento de Córdoba. Debido a los diferentes servicios ecosistémicos que ofrece esta ciénaga en el territorio y a la diversidad que se encuentra en ella, es considerada como uno de los ecosistemas estratégicos del departamento. Sin embargo, a causa de las diferentes actividades antrópicas que se han presentado alrededor de la ciénaga de Betancí se ha presentado una transformación del uso del suelo en la zona como consecuencia de la expansión de la frontera agrícola y ganadera, lo cual ha generado la pérdida casi total del bosque secundario ( Díaz Carvajal , 2017).\par
Por otra parte, Colombia es uno de los principales países latinoamericanos que se caracterizan por sus abundantes recursos naturales presentes en los diferentes ecosistemas del territorio. Los servicios ambientales proveídos por estos ecosistemas han contribuido al desarrollo humano y social de la población colombiana. Sin embargo, las diferentes actividades antrópicas que se ha desarrollado dentro de estos ecosistemas estratégicos han generado una transformación en la cobertura vegetal y el uso del suelo en diferentes zonas, debido a un aprovechamiento inadecuado de los recursos. Este esquema de desarrollo insostenible ha producido grandes impactos ambientales que ponen en riesgo la estabilidad y el bienestar de las presentes y futuras generaciones \hyperref[b20]{(Sánchez Pérez, 2018)}. London Journal of Research in Science: Natural and Formal Los páramos son uno de los ecosistemas estratégicos que presentan un bajo nivel de freático por lo que son considerados muy quebradizos , esto debido a que aproximadamente un 75\% del terreno corresponde a agua (Cuervo Barahona, Cely Reyes, \& Moreno Pérez, 2016). En este contexto, las actividades productivas como la ganadería presentan efectos negativos en la estructura del páramo debido a la pérdida de las propiedades retenedoras de agua y carbono en la matriz del suelo y la compactación de este como producto del pisoteo del ganado (Rivera Ospina, 2013). Los suelos del páramo se caracterizan por contener tres veces más reservas de CO en comparación con la biomasa presente sobre la tierra, además, contienen el doble de carbono que el que se registra en la atmósfera; por tal motivo, la materia orgánica del suelo en este tipo de ecosistemas es una fuente peligrosa para la generación de dióxido de carbono lo cual contribuye a una de las principales problemáticas ambientales como el calentamiento global \hyperref[b8]{(García Portilla, 2003)}.\par
A partir de la fragilidad de los ecosistemas de páramo dentro de territorio colombiano se han generado diferentes normativas para su protección y preservación, como lo es el artículo 173 de la Ley 1753 del 2015, en el cual se estipuló que en estas zonas no se pueden desarrollar actividades productivas que involucren el uso de recursos naturales no renovables. Con base en lo anterior, la población que se encuentra asentada en este tipo de ecosistemas no puede continuar con sus actividades económicas tradicionales. Por lo anterior, es indispensable formular alternativas laborales para las comunidades campesinas residentes en estos ecosistemas ( Rubiano \hyperref[b18]{Galvis, 2015)}.\par
El páramo del Siscunsí, ubicado al sur oriente del municipio de Sogamoso, es uno de los principales ecosistemas estratégico del departamento de Boyacá ya que presenta una extensión de 5.942 ha (Sánchez Pérez, 2018). El principal afluente del páramo es la quebrada La Iglesia, también conocida como quebrada Las Cañas, la cual es uno de los principales drenajes que alimentan al Rio Cusiana dentro del departamento.\par
Considerando lo anteriormente expuesto, es claro que existe la necesidad de conocer el potencial actual de generación de escorrentía en este ecosistema estratégico, y de este modo que las autoridades competentes cuenten con acceso a herramientas e información científica para tener la capacidad de formular acciones de remediación al impacto ambiental provocado sobre este ecosistema de páramo. De esta manera, los objetivos del presente estudio fueron determinar el CN en la microcuenca de la quebrada Las Cañas del páramo de Siscunsí y así reconocer la dinámica hidrológica de este ecosistema en cuanto al potencial de generación de escorrentía, así como identificar las coberturas del suelo a partir de un mosaico multitemporal de imágenes LANDSAT, asignar los grupos hidrológicos del suelo y comparar los resultados de los caudales obtenidos para los diferentes años de estudio a partir del método racional empleando los modelos de Frederich Bell y Dick Peschke para el cálculo de las intensidades, con los resultados entregados por las medias aritméticas analizados a partir del software QGIS de las imágenes ráster de la escorrentía superficial. 
\section[{II. MATERIALES Y MÉTODOS}]{II. MATERIALES Y MÉTODOS} 
\section[{Identificación del cambio de coberturas vegetales y uso del suelo}]{Identificación del cambio de coberturas vegetales y uso del suelo} 
\section[{Delimitación de la cuenca de estudio}]{Delimitación de la cuenca de estudio}\par
Para obtener las imágenes satelitales se realizó la descarga del Modelo de Elevación Digital (DEM, por sus siglas en inglés), obtenido del Servicio de Satélites de Alaska (Alaska Satellite Facility), la cual tiene una resolución espacial por píxel de 12.5 metros. Por medio de la red de drenajes obtenida a partir del procesamiento del DEM y de la información suministrada por el Instituto Geográfico Agustín Codazzi (IGAC), se realizó la delimitación la microcuenca de la quebrada Las Cañas. Estos datos fueron procesados por medio del software QGIS 2.18, donde se realizó la corrección de píxeles, así como la London Journal of Research in Science: Natural and Formal determinación de la dirección y acumulación de flujo y el respectivo análisis de cuenca. En la figura 1, se observa el mapa de localización y delimitación del área de estudio.\par
Figura 1: Localización y delimitación área de estudio 
\section[{Clasificación e identificación del cambio de coberturas y uso del suelo}]{Clasificación e identificación del cambio de coberturas y uso del suelo}\par
Posteriormente, a partir del software QGIS se realizó la identificación del uso del suelo teniendo en cuenta las coberturas vegetales, empleando modelos de predicción. Para la caracterización de uso del suelo se aplicó la metodología Corine Land Cover escala 1:25000 y la información suministrada por el documento Leyenda Nacional de Coberturas de la Tierra Metodología Corine Land Cover escala 1:100000, adaptadas cada una a Colombia. Los datos fueron validados y rectificados con las visitas realizadas a campo.\par
En la figura 2, se presenta el diagrama metodológico para la clasificación e identificación del cambio de coberturas vegetales y uso del suelo.     
\section[{London}]{London} 
\section[{Fuente: Autores}]{Fuente: Autores}\par
Es importante resaltar que se realizó una corrección y normalización de las imágenes satelitales debido a que estas presentaban contrastes diferentes, los cuales pudieron ser producidos por la posición en que se encontraba el sol al momento de su captura y por el relieve de la zona. La corrección topográfica se llevó a cabo por medio de la herramienta "Topographic correction" del software QGIS aplicando los métodos Coseno corrector \hyperref[b3]{(Civco, 1989)}   Posteriormente, se ajustó la clasificación de coberturas terrestres de manera supervisada. Esta clasificación es pertinente por la relación de las clases de cobertura existentes en el área de estudio. La reclasificación de las imágenes satelitales se realizó con la metodología Corine Land Cover adaptada para Colombia mencionada anteriormente \hyperref[b12]{(Menéndez \& Núñez, 2009)}. En este proceso se tuvieron en cuenta, cuatro tipos de cobertura vegetal que fueron identificados en la microcuenca, presentados en la Tabla 1. Estas coberturas se relacionan con los usos de suelo presentes en la región, por medio del POT del municipio de Sogamoso, el EOT del municipio de Aquitania y el POMCA de la cuenca del río Cusiana.\par
Para verificar y corregir la información obtenida se hizo uso en el procesamiento geográfica en QGIS y de la información teórica descrita anteriormente, se realizó una visita al páramo Fuente: Autores identificó cada una de las clases de cobertura y uso del suelo, permitiendo ajustar los datos a las 4 clases de presentadas en la Tabla 1.\par
Tabla 1: Clasificación de cobertura vegetal y uso de suelo  
\section[{Fuente: Autores}]{Fuente: Autores}\par
Con la información obtenida se elaboraron los mapas de uso de suelo de los años 1989,1997, 2002, 2008, 2014 y 2017, por medio de los cuales se llevó a cabo el análisis multitemporal. Se empleó el plugin Molusce del software QGIS para realizar las predicciones de los escenarios de uso del suelo, el cual permite analizar, modelar y simular las transformaciones de la cobertura vegetal y uso del suelo en una zona específica (Villalobos Chacón , 2019).\par
Para la aplicación de esta herramienta es necesario contar con las capas reclasificadas de dos años, así como una combinación de variables que están relacionas a la humedad, erosión y deslizamientos. Estas variables son la distancia a las vías empleando los datos de la red vial (figura 5), drenajes obtenidos del IGAC (figura 6), y también las pendientes dadas en porcentaje que presenta el terreno (figura 7).  En la figura 7, se observa el mapa de pendientes en grados obtenido en QGIS, con la herramienta Slope. 
\section[{Fuente: Autores}]{Fuente: Autores}\par
Figura 7: Mapa pendiente de la cuenca expresada en grados A partir de modelos matemáticos preestablecidos, la herramienta Molusce genera mapas de predicción de los años presentados en la tabla 2. Los periodos obtenidos por medio de los pares dan como resultado una diferencia temporal, la cual, al ser añadida al último año de los pares indica el año a proyectar.\par
Tabla 2: Predicciones realizadas de uso de suelos  
\section[{Fuente: Autores}]{Fuente: Autores}\par
La precisión del modelo se indica en el porcentaje del índice de Kappa. Este es un índice el cual mide la relación entre la clasificación y los píxeles de realidad en el terreno. Los valores de esta constante k indican la concordancia exacta, mientras que un valor de cero no representa ningún tipo de relación (Villalobos Chacón , 2019). 
\section[{Análisis del cálculo del número de curva (CN)}]{Análisis del cálculo del número de curva (CN)}\par
La estimación de la precipitación efectiva se realizó por medio del método propuesto por el Soil Conservation Service (SCS), el cual emplea variables como la cobertura y uso de suelo, las condiciones de humedad antecedente (AMC) y la lluvia acumulada en el área de estudio con el fin de fijar un número de curva (CN) que representa tal interrelación. Para la formulación del CN, se tiene en cuenta que la escorrentía directa (Pe) siempre es menor o igual al agua que precipita en la cuenca (P); de igual manera, después de que comienza la precipitación efectiva, la profundidad adicional del agua retenida (Fa) es menor o igual a la retención potencial máxima (S) (Ven Te, 1984).\par
Con el método de la SCS la escorrentía directa (Q), o precipitación efectiva, es calculada a partir de la siguiente ecuación:\par
(1) ? = (?-0.2?) (?+0.8?) 
\section[{Donde}]{Donde}\par
Q es la escorrentía directa o precipitación efectiva, en pulgadas. P es la precipitación considerada, en pulgadas.\par
S es la diferencia potencial máxima entre P y Q en el momento de inicio de la tormenta y representa proporcionalmente la pérdida de escorrentía por infiltración, intercepción y almacenamiento superficial.\par
Para el cálculo de la diferencia potencial máxima en pulgadas, se aplica la siguiente ecuación.\par
(2) ? = 100 ?? -10 Donde CN es el número de curva. 
\section[{Parámetros del modelo (AMC y Tipo de Suelo Hidrológico)}]{Parámetros del modelo (AMC y Tipo de Suelo Hidrológico)}\par
Las condiciones de humedad antecedente (Antecedent Moisture Conditions, AMC), se agrupan en tres condiciones básicas, como se observa en la Tabla 3. Por otra parte, el tipo de suelo hidrológico se determinó de acuerdo a las condiciones de infiltración subsuperficial que influyen en la precipitación efectiva en la cuenca, teniendo en cuenta el tipo agrológico que se encuentra en la misma. Además, se considera la cobertura vegetal, ya que la vegetación tiene una fuerte incidencia en el agua que escurre dentro de una zona, no solo por establecer el coeficiente de escorrentía sino también por la capacidad que presentan las plantas en absorber el agua en escenarios de humedad.\par
Para la clasificación de suelos se tuvo en cuenta las unidades geológicas del área de estudio identificadas a partir de la información suministrada por la plancha 192 del Servicio Geológico Colombiano, el POT del municipio de Sogamoso y EOT de Aquitania. En la figura 8 se observan las unidades geológicas de la microcuenca de la quebrada las Cañas.   
\section[{Cálculo de caudales a partir de intensidades obtenidas por los modelos de Dick Peschke y Frederich Bell}]{Cálculo de caudales a partir de intensidades obtenidas por los modelos de Dick Peschke y Frederich Bell}\par
Se descargaron los datos pluviométricos de estaciones climáticas del IDEAM, de las cuales 5 son de clase meteorológica (MET) y 6 de clase hidrológica (HID), tal como se observa en la tabla 5. En la figura 10, se representan los polígonos de Thiessen obtenidos en el programa QGIS, donde se observa que la estación que presenta mayor aferencia sobre la microcuenca de estudio es la estación Puente Las Cañas.   Las curvas de intensidad, duración y frecuencia (IDF) se construyeron a partir de los valores de precipitación máxima en 24 horas de las estaciones Toquilla y Puente Las Cañas, durante el periodo de 1971 a 2018. Se llevó a cabo un muestreo estadístico de los datos en intervalos de tiempo de manera discreta y estable, con el objetivo de seleccionar para cada año los datos extremos de precipitación para diferentes valores de tiempo. Posteriormente se llevó a cabo un análisis de frecuencia a través de diferentes tipos de distribuciones tales como Gumbel y log-Gumbel, log-normal de 2 y 3 parámetros, probabilidad normal, gamma de 2 y 3 parámetros, log-Pearson tipo III, y finalmente para calcular los periodos de retorno de 5, 10, 20 y 50 años se realizó una prueba de ajuste de probabilidad Smirnov-Kolmogorov empleando los softwares Excel e Hidroesta 2 .\par
Con los datos de Pmax en 24 horas para los distintos periodos de retorno, se procedió a distribuir las precipitaciones en tiempos inferiores a 60 minutos y de 1 a 24 horas, y de este modo generar las curvas IDF, empleando los siguientes modelos: 
\section[{Modelo de Dick Peschke}]{Modelo de Dick Peschke}\par
Este modelo es usado para el caso de duraciones de tormenta inferiores a 60 minutos, con datos pluviográficos con los que se pueda obtener las intensidades máximas o cuando no se cuente con registros pluviográficos que permitan obtener estas intensidades (Paredes Pinto, 2015). La ecuación que describe este modelo es la siguiente: London Journal of Research in Science: Natural and Formal La intensidad es calcula a partir del cociente entre la precipitación P d y la duración. 
\section[{Modelo de Frederich Bell}]{Modelo de Frederich Bell}\par
De acuerdo con (Paredes Pinto, 2015), el modelo de Frederich Bell es empleado para calcular la lluvia máxima teniendo en cuenta la duración de la tormenta (min), la precipitación máxima de 1 hora de duración y el periodo de retorno de 10 años . La ecuación que describe este modelo es la siguiente: cuenta los valores de precipitación máxima durante 24 h para un periodo de retorno de 10 años, como se observa en la ecuación de Dick Peschke.\par
Por otra parte, la intensidad fue calculada empleando el modelo 2.3.1. 
\section[{Método racional}]{Método racional}\par
Para la determinación de los caudales aplicando el método racional se empleó el Tc de la cuenca (Tiempo de concentración) y la intensidad respectiva a los periodos de retorno (TR) establecidos por medio del método de Dick Peschke y Frederich Bell para estación Puente las Cañas debido a que es la estación que presenta mayor aferencia en el área de estudio (Figura 10). A partir del tiempo de concentración de la microcuenca se realizó la interpolación de la intensidad para los TR correspondientes a 5, 10, 20 y 50 años.\par
La Intensidad de la Lluvia de Diseño se calculó por medio de la Ecuación 6. Es importante resaltar que se consideró la frecuencia como el periodo tiempo de retorno seleccionado y la duración igual al tiempo de concentración de la cuenca ? = Posteriormente, se estableció el coeficiente promedio de escorrentía para el año 2017 por medio del procesamiento del tipo de suelo (Figuras 8 y 9), las pendientes del terreno, la cobertura vegetal obtenida de \hyperref[b11]{(Lemus \& Navarro, 2003)}, los valores geográficos y el área de cada cobertura en el software QGIS. Teniendo en cuenta estas variables se determinó el caudal máximo, como se observa en la siguiente ecuación. 
\section[{? =}]{? =}\par
??? 360 \hyperref[b7]{(7)} Donde Q Es el caudal máximo, en m 3 /s. C Es el coeficiente de escorrentía. I Es la intensidad de lluvia con duración igual al periodo de retorno (mm/h) A Es el área de la cuenca en hectáreas. 
\section[{Comparación de caudales}]{Comparación de caudales}\par
Se llevó a cabo la comparación de los valores de los caudales obtenidos por el método racional con los datos hallados a partir del cálculo de la escorrentía superficial en el software QGIS.\par
Se utilizó el método de estimación inicial del umbral de escorrentía P0 (mm) puesto que no se cuenta con la información suficiente para establecer la lluvia neta; por esta razón, el coeficiente se calculó a partir de las características de la zona de estudio. Por medio del valor inicial de P0 (Tabla 6) se determinó el umbral de escorrentía de acuerdo con el tipo de cobertura vegetal, uso de suelo y los valores de pendientes de la cuenca (Figura 11).  \par
El mapa resultante del procesamiento de estas variables mediante el software SIG se reclasificó de acuerdo con los valores de CN para condiciones normales. Luego, se obtuvieron los mapas de la diferencia potencial máxima por medio del uso de la herramienta "Algebra de mapas". No obstante, a pesar de que en esta metodología emplea el valor de precipitaciones de corta duración para poder evaluar escenarios máximos, en esta investigación se empleó la precipitación total mensual promedio de cada año obtenido a partir de la interpolación realizada con los datos de precipitación registrados por diversas estaciones climáticas en la zona.\par
A partir de los rásteres de la precipitación, CN y factor S se determinó la escorrentía mensual de la cuenca para los años estudiados por medio de la aplicación de la Ecuación 1. Además, mediante los valores promedios de la precipitación efectiva se calculó el volumen de agua. Por último, se determinó el caudal de la cuenca para cada periodo analizado teniendo en cuenta los valores de escorrentía directa y el área de la microcuenca. 
\section[{III. RESULTADOS}]{III. RESULTADOS} 
\section[{Cambio de cobertura y uso del suelo}]{Cambio de cobertura y uso del suelo}\par
Con base a la visita realizada al páramo de Siscunsí, se pudo verificar y correlacionar las 4 categorías de acuerdo con la metodología Corine Land Cover adaptada para Colombia, esto se llevó a cabo con la London Journal of Research in Science: Natural and Formal reclasificación de 15 puntos de observación distribuidos espacialmente a lo largo de la microcuenca Las Cañas. En la Figura 12 se puede evidenciar los puntos de control en la imagen satelital junto con la imagen ajustada y reclasificada en las 4 clases de uso de suelo para el año 2017.\par
Figura 12: Puntos de control en campo junto con imagen en falso color y uso del suelo reclasificado, año 201\par
Por otra parte, a partir del procesamiento de las imágenes satelitales en el software QGIS se pudo establecer los cambios de cobertura de los 4 tipos de uso de suelo que se presentan en el área de estudio durante el periodo 1989-2017. Como se observa en la Tabla 7, las actividades agropecuarias y las áreas para la conservación fueron las categorías que tuvieron un aumento en su extensión dentro de la microcuenca de 49,5 ha y 260,01 ha respectivamente. El incremento de cobertura de las actividades agropecuarias es la consecuencia de la expansión de la frontera agrícola y ganadera durante los 28 años analizados, dando como consecuencia la transformación de otros tipos de usos como los sistemas forestales, los cuales registraron una pérdida de 306,45 ha, en este tipo de uso de suelo. Por otra parte, la ganancia de cobertura de las áreas de conservación es el resultado de las diferentes actividades que se han planteado y desarrollado en la zona a lo largo de tiempo con miras a la conservación del páramo de Siscunsí.\par
En la Figura 13 se ilustran los cambios que presentó la microcuenca Las Cañas en los años 1989-2017 con respecto a las 4 categorías identificadas de uso de suelo. Figura 13: Cambios de uso de suelo durante los años 1989 (Derecha) y 2017 (Izquierda)\par
Por otro lado, empleando las imágenes reclasificadas se realizó la simulación de uso de suelo para los años 2014, 2020, 2026 y 2027 desarrollando la metodología descrita en el apartado anterior. 
\section[{Predicción uso de suelo 2014 a partir de imágenes 2002 y 2008}]{Predicción uso de suelo 2014 a partir de imágenes 2002 y 2008}\par
La proyección de la cobertura de uso de suelo del año 2014 se llevó a cabo con el fin de realizar una comparación de los resultados obtenidos de la simulación (Figura de la izquierda) con la imagen Landsat procesada en el software de este año (Figura de la derecha). Como se observa en la Tabla 8, las actividades agropecuarias y los sistemas forestales son las categorías de uso de suelo que presentaron una diferencia notable teniendo la información obtenida a partir de la simulación y el procesamiento de la imagen satelital. Los valores del área de cobertura dentro de la cuenca de estos dos usos difieren 122 ha y 90 ha respectivamente.\par
Es importante resaltar que la predicción realizada por medio de los resultados del cambio de uso del suelo entre el 2002 y 2008 no tuvo en cuenta la tasa de reforestación que se ha registrado en el páramo, es por esto por lo que los sistemas forestales y las áreas para la conservación son los usos de suelo que presentaron una mayor cobertura en la imagen simulada que en la imagen procesada de la zona de estudio.  El índice de exactitud y el índice Kappa entre estas dos imágenes fue del 73.5\% y 0.624 indicando un valor de concordancia bueno, lo cual significa que el comportamiento calculado y procesado del uso de suelo dentro de la cuenca presenta semejanza por esta razón los modelos de predicción de otros periodos presentaran resultados con una certeza aceptable. presentan crecimiento de sus áreas. Los porcentajes de cambios ilustrados en la Tabla 9 permiten observar que las diferencias de las coberturas totales entre estos dos modelos no son mayores a 3.2\%.\par
Además, al validar estas dos predicciones realizadas para el año 2020 se obtuvo un porcentaje de exactitud de 75.9\% y un índice Kappa del 0.612 reflejando una correlación apropiada entre los datos. Estos valores ratifican que las simulaciones realizadas para los años posteriores aplicando el modelo de esta investigación generaran resultados confiables respecto al uso de suelo en la cuenca.   Los modelos de predicción del uso del suelo para los años 2026 y 2027 difieren significativamente en algunas categorías debido a los distintos periodos empleados para la predicción de estos dos años. Como se observa en la Figura 26, las áreas de conservación en el año 2026 presentaron una mayor área London Journal of Research in Science: Natural and Formal (484,11 ha) en relación con el año 2027 (382,88 ha) debido a que durante los años 2002 y 2014 (Años empleados para la predicción del modelo del 2026) se presentó una mayor expansión en este tipo de suelo dentro de la cuenca que entre los años 1989 y 2008 (Años empleados para la predicción del modelo del 2027). Por otra parte, los sistemas forestales mostraron una tendencia contraria al registrarse una mayor cobertura en el año 2027 (536,03 ha) en comparación con el año 2026 (426,52). En cuanto a los cuerpos de agua y las actividades agropecuarias las diferencias de coberturas de estos usos de suelo dentro de la cuenca en los dos años simulados tan solo difirieron en 0,63 ha y 8,91 ha respectivamente como se observa en la Tabla 10.   Para la verificación , correlación y validación de estos modelos procesados se obtuvo un índice Kappa de 0.46 que representa una concordancia moderada en su desarrollo y una exactitud de 62.8 \% de validación lo que indica que existe similitud en los mapas obtenidos, esto quiere decir, que si el proceso natural se mantiene se tendrán sistemas forestales y áreas de conservación como lo muestra cada mapa London Journal of Research in Science: Natural and Formal de predicción para los años (2026-2027). Asimismo, los índices obtenidos de los modelos 2014 y 2020 arrojaron una concordancia buena, lo cual se debe a que la validación de estos datos se realizó con la imagen satelital del año 2002, lo que genera que el modelo de predicción del año 2026 se asemeje a la realidad. 
\section[{Oferta hídrica de la cuenca empleando la precipitación efectiva}]{Oferta hídrica de la cuenca empleando la precipitación efectiva}\par
De acuerdo con el procesamiento de los datos en el software QGIS, en la Figura 17 se observan los mapas del número de curva obtenidos durante el periodo de estudio (1989-2017) y las predicciones (2020,2026 y 2027).     Teniendo en cuenta los valores promedio de cada una de las imágenes ráster de la precipitación efectiva determinados para el periodo analizado, en la Tabla 11 se presenta el volumen de la microcuenca Las Cañas. 
\section[{London}]{London}\par
Tabla 11: Volumen de agua calculado a partir de los valores promedio de la Precipitación Efectiva  
\section[{Fuente: Autores}]{Fuente: Autores}\par
La oferta hídrica de la microcuenca se calculó por medio de la conversión de los milímetros de escorrentía superficial a las unidades correspondientes de caudal, esto se obtuvo empleando el área total (1799 ha) y los valores máximos de precipitación efectiva para los distintos años estudiados, los resultados se ilustran en la Tabla 12. Mediante el programa Hidroesta2 se estimaron los valores de caudal; al mismo tiempo, por medio del software QGIS se determinó el tipo de uso de suelo, que, junto con los valores de área y precipitaciones de la microcuenca, se estimó el volumen de agua de esta.  
\section[{Comparación de la oferta hídrica empleando la precipitación efectiva y el método racional para el año 2027}]{Comparación de la oferta hídrica empleando la precipitación efectiva y el método racional para el año 2027}\par
Por medio de la Tabla 12 se puede evidenciar que teniendo en cuenta los cambios de uso del suelo y el cálculo de la precipitación efectiva, los caudales obtenidos para los años estudiados no presentan una tendencia uniforme debido a que se registran periodos en los cuales estos valores aumentan mientras que en otros la oferta hídrica de la cuenca disminuye.\par
Por otro lado, los caudales obtenidos por medio del método racional muestran que los resultados obtenidos aplicando el modelo de Frederich Bell aumentan a medida del transcurso de los periodos de retorno analizados, mientras que los valores obtenidos por medio del modelo de Dick Peschke descienden a medida del tiempo (Tabla 16).\par
En cuanto a la comparación de la oferta hídrica para el año 2027, en la Tabla 17 se presentan los resultados obtenidos tanto en la simulación teniendo en cuenta la precipitación efectiva de la microcuenca y los caudales calculados a partir del método racional en el tiempo de retorno de 10 años.\par
Es necesario subrayar que la relación de caudal para el periodo de retorno de 10 años, que para el estudio es el año 2027, con el caudal obtenido por la simulación es de 1:2 por medio del modelo de Frederich Bell y aproximadamente de 1:2.5 empleando el método de Dick Peschke. Para finalizar, teniendo los valores registrados en la Tabla 12, los cuales representan los caudales calculados por medio de la precipitación efectiva, presentan un promedio de 0.46 m 3 /s. No obstante, el promedio de los caudales teniendo en cuenta los datos obtenidos con el método racional empleando el modelo Frederich Bell y Dick Peschke son de 0.247m 3  /s y de 0.  
\section[{IV. CONCLUSIONES}]{IV. CONCLUSIONES}\par
El incremento de áreas en las coberturas vegetales se ha dado significativamente en la microcuenca según los modelos analizados esto en áreas nativas y especies endémicas como lo son los frailejones (Espeletia grandiflora) representando un valor positivo agregado a la investigación y a la zona de estudio donde se puede reconocer el esfuerzo por la protección de estos parques naturales.\par
A partir del análisis multitemporal realizado por medio de software SIG, se logran identificar las transformaciones del uso del suelo durante un periodo de tiempo determinado, además de proveer información actualizada de las zonas en riesgo o con cierto grado de vulnerabilidad, lo cual permite contar con insumos útiles para la toma de decisiones por parte de las autoridades competentes.\par
La variación de los usos de suelos entre los años de 1989 y 2017 indican que los cuerpos de agua han disminuido en un 30\% de su área, así como las actividades agropecuarias que se han visto reducidas en un 11\% del total del área de estudio. Los sistemas forestales han aumentado en un 33\%, mientras que las áreas para la conservación incrementaron en un 30\%, esto gracias a la reforestación y a la población cercana al páramo .\par
Para estimar la oferta hídrica sobre la cobertura y uso de suelo a través del tiempo, se llevó a cabo un análisis multitemporal de estos usos, donde se registran las ganancias y pérdidas de cada uno. El número de curva se evaluó por medio de los caudales máximos empleando los registros históricos de la estación pluviométrica Toquilla y la estación Puente Las Cañas, con un análisis de frecuencia realizado por medio de diferentes distribuciones de probabilidad. La prueba de Kolmogórov-Smirnov reportó que en la estación Puente Las Cañas, la distribución más apropiada fue log-Gumbel (0.1195), mientras que para la estación Toquilla, la distribución normal fue la mejor puesto que presentó un ajuste con un valor de 0.0871.\par
Los resultados de la oferta hídrica empleando la precipitación efectiva y el método racional evidencian que con el método de la SCS el promedio de los valores de caudal fue mayor (0.46 m 3 /s) en comparación con los resultados obtenidos a partir del método racional empleando los modelos de Frederich Bell y Dick Peschke donde se presentaron valores de caudal de 0.247m 3  /s y de 0.162m 3 /s respectivamente. Además, teniendo en cuenta el caudal simulado para un periodo de retorno de 10 años (2027) se identificó que la relación de la oferta hídrica por el método de escorrentía superficial con el método de Frederich Bell es de 1:2 mientras que con el método de Dick Peschke es de 1:2.5.\begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-2.png}
\caption{\label{fig_0}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-3.png}
\caption{\label{fig_1}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{2}\includegraphics[]{image-4.png}
\caption{\label{fig_2}Figura 2 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{3}\includegraphics[]{image-5.png}
\caption{\label{fig_3}Figura 3 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-6.png}
\caption{\label{fig_4}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{4}\includegraphics[]{image-7.png}
\caption{\label{fig_5}Figura 4 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-8.png}
\caption{\label{fig_6}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{5}\includegraphics[]{image-9.png}
\caption{\label{fig_7}Figura 5 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{6}\includegraphics[]{image-10.png}
\caption{\label{fig_8}Figura 6 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-11.png}
\caption{\label{fig_9}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{84}\includegraphics[]{image-12.png}
\caption{\label{fig_10}Figura 8 :Tabla 4 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{9}\includegraphics[]{image-13.png}
\caption{\label{fig_11}Figura 9 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{23}\includegraphics[]{image-14.png}
\caption{\label{fig_12}Fuente: Autores © Volume 23 |}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{10}\includegraphics[]{image-15.png}
\caption{\label{fig_13}Figura 10 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-16.png}
\caption{\label{fig_14}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{10}\includegraphics[]{image-17.png}
\caption{\label{fig_15}10 60}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-18.png}
\caption{\label{fig_16}I}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{6}\includegraphics[]{image-19.png}
\caption{\label{fig_17}Tabla 6 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{11}\includegraphics[]{image-20.png}
\caption{\label{fig_18}Figura 11 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{14}\includegraphics[]{image-21.png}
\caption{\label{fig_19}Figura 14 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-22.png}
\caption{\label{fig_20}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{9}\includegraphics[]{image-23.png}
\caption{\label{fig_21}Tabla 9 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{15}\includegraphics[]{image-24.png}
\caption{\label{fig_22}Figura 15 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{10}\includegraphics[]{image-25.png}
\caption{\label{fig_23}Tabla 10 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{16}\includegraphics[]{image-26.png}
\caption{\label{fig_24}Figura 16 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-27.png}
\caption{\label{fig_25}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{17}\includegraphics[]{image-28.png}
\caption{\label{fig_26}Figura 17 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{18}\includegraphics[]{image-29.png}
\caption{\label{fig_27}Figura 18 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{19}\includegraphics[]{image-30.png}
\caption{\label{fig_28}Figura 19 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{20}\includegraphics[]{image-31.png}
\caption{\label{fig_29}Figura 20 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-32.png}
\caption{\label{fig_30}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-33.png}
\caption{\label{figure33}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-34.png}
\caption{\label{figure34}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{3} \par 
\begin{longtable}{P{0.04175438596491228\textwidth}P{0.3936842105263158\textwidth}P{0.0775438596491228\textwidth}P{0.22070175438596487\textwidth}P{0.11631578947368422\textwidth}}
\tabcellsep ?? ? ( ) =\tabcellsep \multicolumn{2}{l}{4,2 ??(??) 10-0,058?? (??)}\tabcellsep (3)\\
\tabcellsep \multicolumn{2}{l}{?? ??? ( ) =}\tabcellsep 23 ??(??) 0,13?? (??)\tabcellsep (4)\\
\tabcellsep \multicolumn{3}{l}{Lluvia antecedente total de 5 días (pulgadas)}\\
Grupo AMC\tabcellsep \tabcellsep \tabcellsep Estación activa (de\tabcellsep Tipo de condición\\
\tabcellsep \multicolumn{2}{l}{Estación inactiva (seca)}\tabcellsep \\
\tabcellsep \tabcellsep \tabcellsep crecimiento)\\
I\tabcellsep <0,5\tabcellsep \tabcellsep <1,4\tabcellsep Seca\\
II\tabcellsep 0,5 a 1,1\tabcellsep \tabcellsep 1,4 a 2,1\tabcellsep Normal\\
II\tabcellsep Sobre 1,1\tabcellsep \tabcellsep Sobre 2,1\tabcellsep Húmeda\\
\tabcellsep \multicolumn{3}{l}{Information Systems}\end{longtable} \par
 
\begin{quote}
\par
\par
\par
Fuente:\hyperref[b0]{(Avellaneda Córdova , 2020)} Los números de curva se aplican para condiciones de humedad antecedente normales (AMC II). Las otras dos condiciones se determinan a partir de las siguientes ecuaciones:London Journal of Research in Science: Natural and Formal 15 35 Determination of the Hydrological Parameter of the Curve Number (NC) for the Las Cañas Stream in the Siscunsí Paramo by means of Geographic\end{quote}

\caption{\label{tab_2}Tabla 3 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{} \par 
\begin{longtable}{P{0.32320880245649947\textwidth}P{0.01653019447287615\textwidth}P{0.052200614124872056\textwidth}P{0.050460593654042984\textwidth}P{0.023490276356192427\textwidth}P{0.011745138178096214\textwidth}P{0.007395087001023541\textwidth}P{0.050025588536335715\textwidth}P{0.010440122824974411\textwidth}P{0.07395087001023541\textwidth}P{0.01348515864892528\textwidth}P{0.04350051177072672\textwidth}P{0.06177072671443193\textwidth}P{0.06177072671443193\textwidth}P{0.029145342886386896\textwidth}P{0.018705220061412484\textwidth}P{0.0021750255885363357\textwidth}}
\tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \multicolumn{3}{l}{Intensidad (I) en mm/h}\\
\tabcellsep \tabcellsep \tabcellsep \multicolumn{7}{l}{Modelo de Dick Peschke}\tabcellsep \tabcellsep Modelo de Frederich Bell\\
Duración D\tabcellsep \multicolumn{3}{l}{T = 5 Método}\tabcellsep \multicolumn{2}{l}{T= 10}\tabcellsep \multicolumn{3}{l}{T= 20 Tr (5 años)}\tabcellsep \multicolumn{3}{l}{T=50 Tr (10 años) Duración D Tr (20 años) T = 5}\tabcellsep T=10\tabcellsep T=20 Tr (50 años)\tabcellsep T= 50\\
(minutos)\tabcellsep \tabcellsep años\tabcellsep \tabcellsep \multicolumn{2}{l}{años}\tabcellsep \tabcellsep años\tabcellsep \tabcellsep Años\tabcellsep \multicolumn{2}{l}{(minutos)}\tabcellsep años\tabcellsep años\tabcellsep años\tabcellsep años\\
\multicolumn{4}{l}{Frederich Bell}\tabcellsep \tabcellsep \tabcellsep \tabcellsep 43.3\tabcellsep \tabcellsep \tabcellsep 44.9\tabcellsep 46.3\tabcellsep 48.4\\
\tabcellsep \tabcellsep \tabcellsep \tabcellsep \multicolumn{4}{l}{Intensidad (mm)}\tabcellsep \tabcellsep \tabcellsep \tabcellsep Intensidad (mm)\\
10\tabcellsep \multicolumn{5}{l}{493,62 Dick Peschke 603,23}\tabcellsep \tabcellsep 737,18 30.3\tabcellsep \tabcellsep 960,95\tabcellsep 30.1\tabcellsep 10\tabcellsep 195,94 29.8\tabcellsep 225,7\tabcellsep 259,98 29.5\tabcellsep 313,41\\
20\tabcellsep \tabcellsep 293,51\tabcellsep \tabcellsep \multicolumn{2}{l}{358,68}\tabcellsep \tabcellsep 438,33\tabcellsep \tabcellsep 571,38\tabcellsep \tabcellsep 20\tabcellsep 133,51\tabcellsep 153,78 Fuente: Autores 177,14\tabcellsep 213,54\\
30\tabcellsep \tabcellsep 216,54\tabcellsep \tabcellsep \multicolumn{2}{l}{264,63}\tabcellsep \tabcellsep 323,39\tabcellsep \tabcellsep 421,56\tabcellsep \tabcellsep 30\tabcellsep 106,67\tabcellsep 122,87\tabcellsep 141,53\tabcellsep 170,62\\
\multicolumn{12}{l}{40 Por medio de los valores de Intensidad I y el coeficiente de escorrentía determinado para la zona de 174,52 213,27 260,63 339,75 40 90,97 104,78 120,7 145,5}\\
\multicolumn{12}{l}{50 estudio (0.39), se calculó la oferta hídrica de la cuenca teniendo en cuenta los 2 modelos empleados. En 147,63 180,41 220,47 287,39 50 80,4 92,61 106,67 128,6}\\
\multicolumn{10}{l}{60 la Tabla 16 se registran los resultados obtenidos. 128,76 157,35 192,29 250,66}\tabcellsep \tabcellsep 60\tabcellsep 72,68\tabcellsep 83,72\tabcellsep 96,43\tabcellsep 116,25\\
70\tabcellsep \tabcellsep 114,7\tabcellsep \tabcellsep \multicolumn{2}{l}{140,17}\tabcellsep \tabcellsep 171,3\tabcellsep \tabcellsep 223,29\tabcellsep \tabcellsep 70\tabcellsep 66,74\tabcellsep 76,87\tabcellsep 88,55\tabcellsep 106,74\\
80\tabcellsep \tabcellsep 103,77\tabcellsep \multicolumn{9}{l}{126,81 Tabla 16: Caudales calculados a partir del método racional 154,97 202,01 80 61,98 71,39}\tabcellsep 82,24\tabcellsep 99,14\\
90\tabcellsep \tabcellsep 95\tabcellsep \tabcellsep \multicolumn{2}{l}{116,09}\tabcellsep \tabcellsep 141,87\tabcellsep \tabcellsep 184,93\tabcellsep \tabcellsep 90\tabcellsep 58,07\tabcellsep 66,89\tabcellsep 77,05\tabcellsep 92,88\\
100\tabcellsep \tabcellsep 87,78\tabcellsep \tabcellsep \multicolumn{2}{l}{107,27}\tabcellsep \tabcellsep 131,09\tabcellsep \multicolumn{3}{l}{170,88 Caudal (Q) m 3 /s}\tabcellsep 100\tabcellsep 54,78\tabcellsep 63,1\tabcellsep 72,68\tabcellsep 87,62\\
110\tabcellsep \tabcellsep \multicolumn{4}{l}{81,72 Método 99,87}\tabcellsep \tabcellsep \multicolumn{5}{l}{122,05 Tr (5 años) 159,09 Tr (10 años) 110}\tabcellsep 51,97 Tr (20 años)\tabcellsep 59,86 Tr (50 años) 68,95\tabcellsep 83,12\\
120\tabcellsep \tabcellsep \multicolumn{4}{l}{76,56 Frederich Bell 93,56}\tabcellsep \tabcellsep \multicolumn{2}{l}{114,34 0.235}\tabcellsep 149,04\tabcellsep 0.243\tabcellsep 120\tabcellsep 49,52 0.251\tabcellsep 57,04\tabcellsep 0.262\tabcellsep 65,71\tabcellsep 79,21\\
\tabcellsep \tabcellsep \multicolumn{3}{l}{Dick Peschke}\tabcellsep \tabcellsep \tabcellsep \multicolumn{2}{l}{0.164}\tabcellsep \tabcellsep 0.163\tabcellsep 0.161\tabcellsep 0.160\\
\tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep Fuente: Autores\\
\multicolumn{12}{l}{Tabla 12: Caudal promedio calculado a partir de Escorrentía Superficial Tabla 14: Intensidad D empleando los 2 modelos para cada periodo de retorno considerando el Tc}\\
Año\tabcellsep \tabcellsep \multicolumn{2}{l}{1989}\tabcellsep \tabcellsep \multicolumn{2}{l}{1997}\tabcellsep \multicolumn{2}{l}{2002}\tabcellsep \multicolumn{3}{l}{2008 Intensidad (D) en mm/h 2014 2017}\tabcellsep 2020 2026\tabcellsep 2027\\
\multicolumn{2}{l}{Q(m 3 /s)}\tabcellsep \multicolumn{3}{l}{0,521 Método}\tabcellsep \multicolumn{4}{l}{0,486 Tr (5 años) 0,451}\tabcellsep \multicolumn{3}{l}{0,500 Tr (10 años) 0,430}\tabcellsep 0,444 Tr (20 años)\tabcellsep 0,451 Tr (50 años) 0,437\tabcellsep 0,423\\
\tabcellsep \tabcellsep \multicolumn{3}{l}{Frederich Bell}\tabcellsep \tabcellsep \tabcellsep 68.4\tabcellsep \tabcellsep \tabcellsep 78.7\tabcellsep 90.8\tabcellsep Fuente: Autores 109.4\\
\multicolumn{12}{l}{3.3 Oferta hídrica de la cuenca de estudio empleando el método racional Dick Peschke 118.63 144.9 177.2 230.9}\\
\tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep Fuente: Autores\end{longtable} \par
 
\begin{quote}
\par
\par
\par
En la Tabla 13 y Figura 21 se presentan los resultados de intensidad, duración y frecuencia de la estación Puente Las Cañas durante los 4 periodos de retorno (5, 10,20 y 50 años), aplicando los modelos de Dick Peschke y Frederich Bell. Cabe la pena resaltar que estos cálculos también se realizaron para la estación Toquilla, sin embargo, solo se presentan los resultados de Puente Las Cañas por ser la que mejor representa las condiciones de la cuenca.London Journal of Research in Science: Natural and Formal 15 51 Determination of the Hydrological Parameter of the Curve Number (NC) for the Las Cañas Stream in the Siscunsí Paramo by means of Geographic Information Systems Tabla 13: Fuente: Autores Fuente: Autores Figura 21: Curvas IDF aplicando los modelos Dick Peschke y Frederich BellCon el fin de determinar los valores de intensidad durante el tiempo de concentración de la cuenca (67,2 minutos), se realizó una interpolación para cada uno de los periodos de retorno empleando los 2 modelos descritos anteriormente. Los resultados se registran en la Tabla 14.En la Tabla 15 se presentan los resultados de la intensidad I calculada por medio de los valores de intensidad D.\end{quote}

\caption{\label{tab_9}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{} \par 
\begin{longtable}{P{0.5862068965517242\textwidth}P{0.24620689655172415\textwidth}P{0.017586206896551722\textwidth}}
Método\tabcellsep Caudal (m 3\tabcellsep /s)\\
Método Racional-Modelo Frederich Bell\tabcellsep 0.243\tabcellsep \\
Método Racional-Modelo Dick Peschke\tabcellsep 0.163\tabcellsep \\
Precipitación Efectiva\tabcellsep 0.423\tabcellsep \\
\tabcellsep \multicolumn{2}{l}{Fuente: Autores.}\end{longtable} \par
 
\caption{\label{tab_10}}\end{figure}
 			\label{foot_0}\footnote{\label{foot_0} Volume 23 | Issue 2 | Compilation 1.0 © 2023 Great Britain Journal Press} 			\label{foot_1}\footnote{\label{foot_1} Volume 23 | Issue | Compilation 1.0 © 2023 Great Britain Journal Press} 		 		\backmatter  			  				\begin{bibitemlist}{1}
\bibitem[ London Journal of Research in Science: Natural and Formal]{b5}\label{b5} 	 		\textit{},  	 	 		\textit{London Journal of Research in Science: Natural and Formal}  		 	 
\bibitem[García Portilla ()]{b8}\label{b8} 	 		\textit{},  		 			J García Portilla 		.  		2003. Bogotá, D.C.  		 			Universidad Javeriana 		 	 
\bibitem[Calvache ()]{b1}\label{b1} 	 		\textit{},  		 			M Calvache 		.  		2010.  	 
\bibitem[Rivera Ospina ()]{b15}\label{b15} 	 		\textit{},  		 			D Rivera Ospina 		.  		2013.  	 
\bibitem[Sánchez Pérez ()]{b20}\label{b20} 	 		\textit{},  		 			D F Sánchez Pérez 		.  		2018.  	 
\bibitem[Romano Castillo ()]{b17}\label{b17} 	 		\textit{Algoritmo computacional para describir la dinámica del ciclo hidrológico superficial en una cuenca a partir de cartas temáticas topográficas, edafológicas, uso de suelo y vegetación},  		 			E A Romano Castillo 		.  		2007. Foro-Red-Mat: Revista electrónica de contenido matemático.  	 
\bibitem[Avellaneda Córdova ()]{b0}\label{b0} 	 		\textit{Análisis comparativo de los métodos Racional modificado Témez, Hidrogramas unitarios SCS, Clark y Snyder en la obtención de caudales máximos para las subcuencas Cañad y Alto Chancay Lambayeq Perú},  		 			A D Avellaneda Córdova 		.  		2020. p. .  		 			Universidad Nacional Pedro Ruiz Gallo 		 	 
\bibitem[Villalobos Chacón ()]{b23}\label{b23} 	 		‘Análisis del potencial de escurrimiento en la microcuenca del río Tiribí entre el 2005 al 2029 para la conservación de los servicios ecosistémicos de regulación hídrica. Heredia: Universidad Nacional de Costa Rica’.  		 			A Villalobos Chacón 		.  	 	 		\textit{London Journal of Research in Science: Natural and Formal}  		2019. 15 p. 55.  	 	 (Determination of the Hydrological Parameter of the Curve Number (NC) for the Las Cañas Stream in the Siscunsí Paramo by means of Geographic Information Systems) 
\bibitem[Cuervo Barahona et al. ()]{b4}\label{b4} 	 		‘Determinación de las fracciones de carbono orgánico en el suelo del páramo la Cortadera’.  		 			E Cuervo Barahona 		,  		 			G Cely Reyes 		,  		 			D Moreno Pérez 		.  	 	 		\textit{Boyacá. Ingenio Magno}  		2016. 7 p. .  	 
\bibitem[Díaz Carvajal ()]{b7}\label{b7} 	 		\textit{Determinación del número de curva en la subcuenca de Betancí (Córdoba, Colombia) mediante teledetección y SIG},  		 			Á Díaz Carvajal 		.  		2017. Ingeniería y Desarrollo. p. .  	 
\bibitem[Castro López ()]{b2}\label{b2} 	 		\textit{El sistema de información geográfica aplicado a la investigación urbana en México: balance y formación técnica},  		 			J Castro López 		.  		2011. México.  		 			Universidad Autónoma Metropolitana Unidad Azcapotzalco 		 	 
\bibitem[Deshmukh et al. ()]{b6}\label{b6} 	 		‘Estimation and comparision of curve numbers based on dynamic land use land cover change, observed rainfall-runoff data and land slope’.  		 			D S Deshmukh 		,  		 			U C Chaube 		,  		 			A Ekube Hailu 		,  		 			D Aberra Gudeta 		,  		 			Tegene Kassa 		.  	 	 		\textit{Journal of Hydrology}  		2013. 492 p. .  	 
\bibitem[Galvis ()]{b18}\label{b18} 	 		 			Rubiano Galvis 		,  		 			S 		.  		\textit{Protección de páramos y derechos campesinos. Aportes jurídicos y de política Subdirección de Servicios Científicos y Proyectos Especiales . Bogotá, D.C: Instituto de Investigación de Recursos Biológicos Alexander von Humboldt},  		2015.  	 
\bibitem[Paredes Pinto ()]{b14}\label{b14} 	 		\textit{Generación de intensidades de precipitación máxima en estaciones sin información pluviográfica (cuenca del río Illpa-Puno)},  		 			R E Paredes Pinto 		.  		2015. Puno-Perú: Universidad Nacional del Altiplano.  	 
\bibitem[Cruz ()]{b21}\label{b21} 	 		\textit{Gestión integrada del riesgo de inundaciones en Colombia},  		 			Sedano Cruz 		,  		 			RK 		.  		2013. 37 p. .  		 			Universitat Politècnica de València 		 	 
\bibitem[Lemus and Navarro ()]{b11}\label{b11} 	 		\textit{Manual para el desarrollo de obras de conservación de suelo},  		 			M Lemus 		,  		 			V G Navarro 		.  		2003. San Fernando (Chile: Corporación Forestal Nacional.  	 
\bibitem[Menéndez and Núñez ()]{b12}\label{b12} 	 		 			M Menéndez 		,  		 			V Núñez 		.  		\textit{El uso de los sensores remotos en los Recursos Naturales Primera parte: La fotografía aérea y la fotointerpretación Argentina},  		2009. 3 p. .  	 
\bibitem[Palomino ()]{b19}\label{b19} 	 		\textit{Modelación precipitación escorrentía en cuencas urbanas-caso de estudio: Distrito de Carmen Alto},  		 			Sánchez Palomino 		.  		2015.  		 			Universidad Nacional de San Cristóbal de Huamanga 		 	 
\bibitem[Ramos ()]{b9}\label{b9} 	 		 			Lao Ramos 		,  		 			B 		.  		\textit{La teledetección y los Sistemas de Información Geográfica para el manejo de las tierras. Revista Ciencias Técnicas Agropecuarias},  		2018. p. .  	 
\bibitem[Rodríguez et al. ()]{b16}\label{b16} 	 		 			A Rodríguez 		,  		 			F Gaspari 		,  		 			G Sinisterra 		,  		 			M Delgado 		,  		 			S Besteiro 		.  		\textit{Evaluación del efecto de la restauración agro-hidrológica mediante la aplicación del modelo hidrológico GeoQ. Revista Investigaciones Agropecuarias},  		2012. p. .  	 
\bibitem[Law and Nichol ()]{b10}\label{b10} 	 		‘Topographic correction for differential illumination effects on IKONOS satellite imagery’.  		 			K H Law 		,  		 			J Nichol 		.  	 	 		\textit{Int. Arch. Photogramm. Remote Sens. Spat. Inform. Sci}  		2004. 35 p. .  	 
\bibitem[Civco ()]{b3}\label{b3} 	 		‘Topographic normalization of Landsat Thematic Mapper digital imagery’.  		 			D L Civco 		.  	 	 		\textit{Photogrammetric Engineering \& Remote Sensing}  		1989. 55  (9)  p. .  	 
\bibitem[Morad ()]{b13}\label{b13} 	 		‘Una aproximación a las ventajas y dificultades de su aplicación’.  		 			M Morad 		.  	 	 		\textit{Sistemas de Información Geográfica y modelizaciones hidrológicas},  		2001. 31 p. .  	 
\bibitem[Ven Te ()]{b22}\label{b22} 	 		 			C Ven Te 		.  		\textit{Handbook of Applied Hydrology},  		1984. McGraw Hill Book Company.  	 
\end{bibitemlist}
 			 		 	 
\end{document}
